UNCLASSIFIED 


_ AD  NUMBER _ 

AD904665 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  SEP  1972.  Other 
requests  shall  be  referred  to  Air  Force  Flight 
Dynamics  Lab.,  Wright-Patterson  AFB,  OH  45433. 


_ AUTHORITY 

AFWAL  ltr  14  Aug  1980 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  SEEN  DELIMITED 
AND  CLEARED  FOR  PUILIC  RELEASE 

UNDER  DOD  DIRECTIVE  5200,20  AND 

NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  POR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


AD904665 


( 


AFFDL-TR-72-40 
Volume  III 


V/STOL  DYNAMICS  AND  AEROELASTIC 
ROTOR-AIRFRAME  TECHNOLOGY 

Volume  III.  User’s  Manuals  —  Computer 
Programs  for  Aeroelastic  Stability  Analysis 
and  Aeroelastic  Prop/Rotor  Loads  Analysis 

A.  K.  Amos 
F.  J.  Tarzanin 

The  Boeing  Company,  Vertol  Division 


TECHNICAL  REPORT  AFFDL-TR-72-40, 
Volume  III 

SEPTEMBER  1972 

4 

Distribution  llmltad  to  U.S.  Gov.rnm.nt  agendas  only; 
test  and  .valuation;  statamant  appllad  18  April  1972. 

Othar  r .quarts  for  this  documant  must  ba  rafarred  to 
tha  AF  Flight  Dynamics  Laboratory,  (FY),  Wright-Pattarson 
AFB,  Ohio  45433. 


Air  Force  Flight  Dynamics  Laboratory 
Air  Force  Systems  Command 
Wright-Patterson  Air  Force  Base,  Ohio 


When  Government  drawings,  specifications,  or  other  data  are 
used  for  any  purpose  other  than  in  connection  with  a  definitely 
related  Government  procurement  operation,  the  United  States 
Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formu¬ 
lated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  impli¬ 
cation  or  otherwise  as  in  any  manner  licensing  the  holder  or 
any  other  person  or  corporation,  or  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 


Copies  of  this  report  should  not  be  returned  unless  return  is 
required  by  security  considerations,  contractual  obligations, 
or  notice  on  a  specific  document. 


AIR  FORCE/56780/2  November  1972-  100 


V/STOL  DYNAMICS  AND  AEROE1ASTIC 
ROTOR-AIRFRAME  TECHNOLOGY 


Volume  III.  User's  Manuals  -  Computer  Programs 

for  Aeroelastic  Stability  Analysis  and 
Aeroelastic  Prop/Rotor  Loads  Analysis 


A.  K.  Amos 
F.  J.  Tarzanin 


The  Boeing  Company,  Vertol  Division 


Distribution  limited  to  U.  S.  Government  agencies  only; 
test  and  evaluation;  statement  applied  18  April  1972. 
Other  requests  for  this  document  must  be  referred  to 
the  AF  Flight  Dynamics  Laboratory,  (FY) ,  Wright  - 
Patterson  Air  Force  Base,  Ohio  45433. 


v  r*j  j 


mL* a 


FOREWORD 

This  report  was  prepared  by  The  Boeing  Company,  Vertol 
Division  of  Philadelphia,  Pennsylvania,  for  the  Aerospace 
Dynamics  Branch,  Vehicle  Dynamics  Division,  Air  Force  Flight 
Dynamics  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio, 
under  Contract  F33615-71-C-1310 .  This  research  is  part  of  a 
continuing  effort  to  develop  new  and  improved  techniques  for 
defining  dynamic  and  aeroelastic  phenomena  for  rotor/propeller- 
powered  V/STOL  flight  vehicles  under  the  Air  Force  Systems 
Command's  exploratory  development  program.  This  contract  was 
initiated  under  Project  1370,  "Dynamic  Problems  in  Military 
Flight  Vehicles,"  Task  137005,  "Prediction  and  Control  of 
Flight  Vehicle  Vibration."  Mr.  A.  R.  Basso  of  the  Aerospace 
Dynamics  Branch  was  the  Project  Engineer. 

The.  final  report  is  presented  in  three  volumes.  The  first 
volume  contains  a  state-of-the-art  review  of  stability  and 
blade  vibratory  loads  in  V/STOL  aircraft.  The  second  volume 
contains  the  development  of  the  analytical  methods,  the  corre¬ 
lation  of  analytical  results  with  experimental  data,  and  the 
results  of  parametric  investigations.  The  third  volume  con¬ 
tains  a  user's  guide  to  the  digital  computer  programs  including 
input  and  output  formats.  The  third  volume  is  not  being 
distributed;  however,  it  is  available  upon  request  from  the 
Air  Force  Flight  Dynamics  Laboratory/FYS ,  Wright-Patterson 
Air  Force  Base,  Ohio  45433. 

Mr.  H.  R.  Alexander  was  The  Boeing  Company,  Vertol  Division 
Project  Engineer. 

This  report  covers  work  conducted  from  February  1971  through 
February  1972.  The  manuscript  was  released  by  the  authors  in 
February  1972  for  publication  as  an  AFFDL  Technical  Report. 

This  Technical  Report  has  been  reviewed  and  is  approved. 

walter(/j.  I^yKYTOW 
Assistant  for  Research  and 
Technology 

Vehicle  Dynamics  Division 
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ABSTRACT 


This  report  provides  user’s  instructions  for  two  computer 
programs,  one  for  aeroelastic  stability  analysis  and  one  for 
aeroelastic  prop/rotor  loads  analysis.  These  programs  are 
commonly  identified  by  Boeing-Vertol  as  C-39  and  C-70,  respec¬ 
tively.  Each  program  is  carefully  described  and  explained, 
symbols  and  sign  conventions  are  identified,  and  input  and 
output  data  are  presented.  A  sample  program  run  for  each 
analysis  is  then  given.  In  addition,  sample  programs  for 
subroutines  D-01  and  A-97,  which  are  used  in  support  of  C-70, 
are  provided.  Notes  and  suggestions  on  program  usage  are 
presented . 


iii 


luVkii 


•PRECEDING  PAGE  BLANK-NOT  FILMED.* 


I 


J 


TABLE  OF  CONTENTS 


p.*g.g. 


Part  I.  User's  Manual  for  Computer  Program  for 

Aeroelastic  Stability  Analysis  .  1 

INTRODUCTION  .  3 

PROGRAM  FEATURES  .  5 

SYSTEM  DEFINITION  .  5 

MODES  AND  SYSTEM  FREEDOMS .  5 

LOOPING  .  6 

INPUT  DATA .  7 

PROGRAM  OUTPUT .  41 

SAMPLE  PROBLEM  .  48 

INPUT  DATA .  4  8 

PROGRAM  OUTPUT . .  .  .  55 

DISCUSSION  OF  PROGRAM  USAGE  .  96 

Part  II.  User's  Manual  for  Computer  Program  for 

Aeroelastic  Prop/Rotor  Loads  Analysis  ....  97 

INTRODUCTION  .  99 

BLADE  PHYSICAL  PROPERTIES .  99 

BLADE  AIRFOIL  CHARACTERISTICS . 100 

TRIM  DATA . 100 

SPECIFICATION  OF  PROGRAM  CONTROLS  .  100 

AIRFOIL  DECKS  .  103 

PROGRAM  INPUT . 120 

SIGN  CONVENTIONS . 150 

DESCRIPTION  OF  OUTPUT  DATA . 155 

SAMPLE  PROGRAM  INPUT  AND  OUTPUT  .  175 

SAMPLE  BLADE  PROPERTY  IDEALIZATION  PROGRAM 

INPUT  AND  OUTPUT . 211 

SAMPLE  TANDEM  HELICOPTER  TRIM  AND  STABILITY 

ANALYSIS  PROGRAM  INPUT  AND  OUTPUT  .  247 


V 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1  Aircraft  Coordinate  Axes  (Fixed-Coordinate 

System) . 43 

2  Schematic  Diagram  of  Rotor  Hub  Showing 

Blade  Mount  Eccentricities  .  44 

3  Rotor  and  Aircraft  Coordinate  Axes  .  45 

4  Local  Coordinate  Axes  for  Typical  Segment 

of  Various  Airframe  Components  .  46 

5  Typical  Landing  Gear  Representation  .  47 

6  Lift  Coefficient  Versus  Angle  of  Attack 

for  All  Airfoils . 105 

7  Section  Drag  Characteristics  at  Large 

Angles  of  Attack . 107 

8  Section  Pitching  Moment  Characteristics 

at  Large  Angles  of  Attack . 109 

9  Coefficient  of  Lift  for  NACA  0011  (Mod) 

Airfoil  Section  .  117 

10  Coefficient  of  Drag  for  NACA  0011  (Mod) 

Airfoil  Section  .  118 

11  Coefficient  of  Pitching  Moment  for  NACA  0011 

(Mod)  Airfoil  Section . .  119 

12  Circulation  Strength  Factor  at  Blade  Tip  .  .  .  137 

13  Circulation  Strength  Factor  at  Blade  Root  .  .  .  138 

14  Fuselage  Geometry  .  150 

15  Vortex  Geometry . 150 

16  Blade  Airfoil  Geometry  .  151 

17  Airfoil  Sign  Convention . 151 

18  Radial  Thrust  Component  Sign  Convention  ....  152 

19  Radial  Force  Component  From  Flap 

Foreshortening  Sign  Convention  .  152 


vi 


Figure  Page 

20  Coriolis  and  Horizontal  Loads  .  152 

21  Coupled  Flap-Pitch  Sign  Convention  .  153 

22  Uncoupled  Lag  Sign  Convention . 153 

23  Rotating  and  Fixed-System  Hub  Loads  Sign 

Convention  (Shaft  Axis  System)  .  154 

24  Steady  Loads  Relative  to  Airstream  (Wind 

Axis  System)  . 3  54 

25  Spanwise  Weight  Distribution  of  CH-47 

Rotor  Blade . 216 

26  Effective  Centroidal  Axis  Location  of  CH-47 

Rotor  Blade . 217 

27  Effective  Polar  Mass  Moment  of  Inertia  About 

Pitch  Axis  of  CH-47  Rotor  Blade . 218 

28  Effective  Flapwise  Stiffness  Distribution  of 

CH-47  Rotor  Blade . 219 

29  Effective  Chordwise  Stiffness  Distribution  of 

CH-47  Rotor  Blade . 229 

30  Torsional  Stiffness  Distribution  of  CH-47 

Rotor  Blade . 221 

31  Effective  Axial  Stiffness  Distribution  of 

CH-47  Rotor  Blade . 222 

32  Effective  Neutral  Axis  Location  of  CH-47 

Rotor  Blade . 223 


vii 


Part  I.  User’s  Manual:  Computer  Program 
for  Aeroelastic  Stability  Analysis 


A.  K.  Amos 


VT  7TT 


r  ..4  - 


F 


vr 


•PP1CEDIM3  PAGE  BLANK-NOT  FUMED. 


INTRODUCTION 

The  program  is  designed  to  generate  and  solve  a  set  of  equations 
representing  flutter  and  vibration  states  of  dynamic  equilibrium 
of  a  general  rotored  aircraft  system  following  an  arbitrary  per¬ 
turbation  from  a  known  steady  state  configuration.  In  matrix 
notation,  these  equations  are  represented  by  . 


[-  [M]  +  ix  (  [C]  +  rz  [AX] )  +  x2([K]  +  rR  [Ar]  )  ]  {q}  =  {o} 

where  [M] ,  [C] ,  and  [K]  are  generalized  inertia,  damping  and 

stiffness  matrices  arising  from  the 
structural  dynamics  of  the  entire  system 
and  the  aerodynamics  of  the  rotor  blades; 

[A  ]  and  [A  1  are  coefficient  matrices  arising  from  the 
R  1  aerodynamics  of  the  airframe  lifting  surfaces; 

{q}  is  a  vector  of  system  freedoms  (generalized  coordin¬ 
ates)  ; 

X  is  a  complex  scalar  parameter; 
rR  and  Tj  are  real  scaling  factors ;  and 
i  is  the  imaginary  number  /-I  . 


The  aircraft  model  on  which  the  program  is  based  comprises  three 
major  subsystems  identified  as  (i)  Rotor,  (ii)  Airframe,  and  (iii) 
Landing  Gear  Subsystems. 

The  rotor  subsystem  is  made  up  of  two  rotor-nacelle  units.  Each 
rotor  has  three  or  more  non-art iculated  flexible  blades  attached 
to  a  teetering  hub.  Steady  state  rotor  configurations  are  pre¬ 
sumed  to  involve  large  blade  deflections  in  flap,  lag  and  twist, 
and  any  arbitrary  inclination  of  the  shaft  axis  with  the  direction 
of  flight.  This  last  feature  permits  treatment  of  tilt  rotor  con¬ 
figurations  ranging  from  hover  to  forward  flight  conditions.  Phy¬ 
sically  the  rotors  may  be  supported  anywhere  in  the  airframe  and 
coordinates  are  provided  for  specifying  the  support  locations  for 
each  application.  Thus  both  helicopter  and  rotored  conventional 
aircraft  systems  are  represented  by  the  general  model  of  the  pro¬ 
gram. 

The  airframe  subsystem  has  the  following  components:  -  right  wing, 
left  wing,  fuselage,  right  tail,  left  tail,  and  fin.  Except  for 
the  fuselage,  all  the  components  are  treated  as  lifting  surfaces. 
The  steady  state  airframe  configuration  is  specified  by  sets  of 
coordinates  capable  of  representing  a  wide  variety  of  practical 
aircraft  systems.  For  example  wing  "tilt"  and  "sweepback"  angles 
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of  aircraft  systems  having  these  features  will  be  reflected  in  the 
coordinates  for  the  wing  components.  Wing  or  tail  dihedral  angle 
has,  however,  not  been  provided  for. 


The  landing  gear  subsystem  consists  of  four  linearized  oleo-strut- 
tire  units.  Their  points  of  attachment  in  the  airframe  are 
identified  by  sets  of  coordinates  to  be  specified  for  the  particu¬ 
lar  systems  being  analyzed. 

The  analytical  basis  of  the  program  lies  in 


(i)  Modal  representation  of  system  deviations  from  steady 
state  equilibrium  configurations; 

(ii)  A  basically  conservative  system  under  such  deviations; 

(iii)  Two-dimensional  quasi-static  aerodynamic  theory  for  the 
rotor  blades; 

(iv)  Two-dimensional  oscillatory  aerodynamic  theory  for  the 
airframe  lifting  bodies;  and 

(v)  Linearly  visco-elastic  oleo  struts  and  tires. 
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PROGRAM  FEATURES 


SYSTEM  DEFINITION 

The  aircraft  model  previously  described  has  been  devised  to  furnish 
a  general  framework  representative  of  a  large  class  of  practical 
systems.  To  this  end  the  system  components  and  their  steady  state 
configurations  have  been  characterized  by  general  coordinates  with 
wide  ranges  of  specification.  In  addition,  a  number  of  bi-valued 
parameters  (option  indices)  have  been  provided  for  eliminating  from 
or  retaining  in  the  model  entire  components  or  subsystems  in  adapt¬ 
ing  it  to  specific  cases. 

Each  of  the  two  rotors  can  be  retained  or  deleted  to  provide  repre¬ 
sentation  of  systems  with  fewer  than  two  rotors.  For  programming 
purposes  the  rotors  have  been  numbered  1  and  2  respectively  with 
the  first  having  a  higher  priority  of  retention  over  the  second. 

Thus  for  a  single  rotored  system  ROTOR  #1  (rather  than  #2)  MUST  be 
specified  as  RETAINED  (see  "Input  Data"  Section) . 

The  airframe  components  can  each  be  specified  as  present  or  not 
present  depending  on  the  particular  system  being  analyzed.  A  half 
aircraft  analysis,  for  example,  is  achieved  by  specifying  one  wing 
and  one  tail  (both  being  either  right  or  left)  as  absent,  and  by 
supplying  fuselage  and  fin  characteristics  representative  of  half 
the  items. 

The  actual  number  (not  exceeding  four)  of  landing  gears  in  the 
particular  system  being  analyzed  is  specified  as  input  to  the  pro¬ 
gram  which  then  limits  the  analysis  to  just  that  many  gears. 

MODES  AND  SYSTEM  FREEDOMS 

Perturbations  of  the  system  about  its  steady  state  equilibrium 
configuration  are  represented  by  a  number  of  modes  falling  into 
three  basic  categories.  The  first  is  made  up  of  BLADE  VIBRATORY 
MODES.  They  are  characterized  by  flap,  lag  and  twist  components 
at  all  points  of  the  blades;  frequencies  of  natural  vibration  at 
specified  rotational  speeds;  and  optionally,  some  measure  of  modal 
damping.  Each  mode  is  associated  with  a  COLLECTIVE  and  IP  CYCLIC 
ROTOR  FREEDOMS.  There  are  thus  three  rotor  freedoms  per  blade  mode. 

AIRFRAME  MODES  constitute  the  second  category.  They  are  defined  for 
each  point  of  the  system  including  the  rotor  blades  (after  blade 
modes  have  been  accounted  for) .  Each  is  associated  with  a  single 
airframe  modal  freedom.  The  third  category  of  modes  arises  from  the 
landing  gears  with  each  mode  defined  to  produce  no  displacements 
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anywhere  except  at  one  of  the  gear  sprung  masses  where  it  gives  rise 
to  a  vertical  displacement  of  unit  magnitude.  Each  of  the  gears  is 
thus  affected  by  only  one  mode  and  hence  one  degree  of  freedom. 

LOOPING 

The  program  features  automatic  cyclirlg  (looping)  through  the  analy¬ 
sis  for  several  problem  cases  representing  parametric  variations  in 
some  of  the  system  characterizations.  In  usage,  a  full  set  of  in¬ 
put  data  is  supplied  for  the  first  case  and  only  updates  of  a  group 
of  pre-established  data  items  for  all  subsequent  cases.  Five  such 
groups  of  data  as  described  below  are  available  for  optional  selec¬ 
tion. 

OPTION  1:  Aircraft  velocity  components  as  well  as  air  density  and 
the  speed  of  sound  vary  from  case  to  case. 

OPTION  2:  Rotor  RPM  and  steady  state  blade  deflections  for  all 
rotors  present  are  updated  for  each  case. 

OPTION  3:  The  items  in  options  1  and  2  together  are  updated  for 
each  case. 

OPTION  4:  Updated  items  consist  of  the  modal  frequencies  and  damp¬ 
ing  factors  for  the  elastic  modes  of  the  system. 

OPTION  5:  The  landing  gear  properties  are  updated  for  each  case. 


INPUT  DATA 


The  program  accepts  data  on  cards  only.  All  data  items  must  be 
based  on  a  slug-foot-second  set  of  units.  However,  rotational 
velocities  and  frequency  items  are  input  as  either  RPM  or  radians 
per  second  as  indicated  later  in  this  section. 


The  input  data  has  been  arranged  in  groups  as  described  below: 

Data  Group  A:  Program  Definition:  Specifies  the  problem  title, 

the  six  components  of  rigid  displacement,  and  a 
series  of  logical  data  to  define  the  general 
features  of  the  current  program  application. 

This  group  is  required  for  all  runs. 


Data  Group  BI:  Rotor  Configuration  Definition:  includes  param¬ 

eters  for  defining  the  rotor  subsystem  as  to  the 
number  of  rotors  involved,  their  mode  of  opera¬ 
tion,  and  the  blade  mode  freedoms  to  be  considered 
in  the  analysis.  Each  of  the  three  modal  freedoms 
(collective,  cyclic  yaw,  cyclic  pitch)  for  any  of 
the  blade  modes  provided  can  be  arbitrarily  elim¬ 
inated  from  the  analysis  where  this  is  desirable 
in  obtaining  a  realistic  representation  of  the 
rotor  system. 

The  group  must  be  specified  for  all  runs. 

Data  Group  BIT:  Rotor  Data  Set  No.  1:  This  group  is  not  required 

if  no  rotors  are  involved  in  the  analysis.  It 
contains  numerical  data  for  all  the  characteris¬ 
tics  of  the  first  (and  for  identical  rotors,  the 
second)  rotor  in  the  system. 

Data  Group  Bill:  Rotor  Data  Set  No.  2:  This  group  of  data  is 

required  only  when  two  dissimilar  rotors  are 
involved  in  the  analysis,  it  is  not  required 
for  a  single  rotor  or  a  pair  of  geometrically 
identical  and  operationally  similar  rotors.  In 
the  later  case  the  first  set  applies  to  both 
rotors. 


Data  Group  BIV:  Nacelle  No.  1  Parameters:  Geometrical  and  basic 

dynamic  properties  of  the  first  nacelle  in  the 
system  are  here  defined.  It  is  presumed  that 
this  nacelle  and  Rotor  No.  1  are  parts  of  the 
same  physical  unit  and  a  nacelle-rotor  coupling 
between  them  is  automatically  generated  using  data 
specified  separately  for  each. 
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Data  Group  BV: 


Nacelle  No.  2  Parameters:  Repetition  of  the 
previous  group  for  the  second  nacelle  if  present. 
Otherwise  this  group  must  be  eliminated. 


Data  Group  Cl: 


Airframe  Configuration  Definition:  Required  for 
all  runs,  it  specifies  the  airframe  components 
present  in  the  analysis. 


Data  Group  CII: 


Right  Wing  Parameters:  Required  if  this  compo¬ 
nent  has  been  specified  present  in  Cl.  Other¬ 
wise  omit  entire  group. 


Data  Groups  CIII 
Through 
CVII: 


Parameters  for  Other  Airframe  Components:  Same 
items  (with  some  modifications)  as  defined  for 
Group  CII  but  applicable  to  the  other  airframe 
components  in  order.  However,  they  are  included 
in  the  data  package  only  if  the  components  they 
apply  to  have  been  specified  in  Cl  as  present. 


Data  Group  CVII I: 


Airframe  Modal  Frequencies  and  Damping  Factors: 
The  modal  frequencies  of  the  airframe  elastic 
modes  and  the  corresponding  damping  factors  due 
to  the  airframe  structure  are  specified  here. 
Entire  group  must  be  eliminated  if  no  elastic 
modes  are  present. 


Data  Group  D: 


Landing  Gear  Properties;  Not  required  if  no 
landing  gears  are  involved  in  the  analysis. 


Data  Group  E: 


Airframe  Aerodynamics  Supplemental  Data:  Com¬ 
prises  sweepback  angles  for  the  airframe  lifting 
bodies  and  a  reduced  frequency  for  the  first 
solution  case. 


When  airframe  aerodynamics  are  being  ignored  in 
the  analysis,  this  group  of  data  must  be  elim¬ 
inated. 


Data  Groups  FI, 
FII. . . FV: 


Data  Updates  for  Looping  Options  1  Through  5: 
The  group  of  data  appropriate  to  the  looping 
option  previously  selected  is  repeated  for  each 
solution  case  following  the  first. 


Following  is  a  tabulation  and  description  of  the  input  data  items 
and  their  governing  input  formats. 
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ITEM 

DESCRIPTION 

UNITS 

4 

Mass  per  unit  length 

slugs/ft 

5 

C.G.  location  aft  of  elastic  axis 

6 

Mass  moment  of  inertia  per  unit  length  about  local 
chord  axis  (assumed  to  pass  through  elastic  center) 

slug-flf/ft 

7 

Mass  moment  of  inertia  per  unit  length  about  normal  to 
chord  axis  passing  through  elastic  center 

slug-ft^/ft 

8 

Average  chord  length  over  segment  length 

ft. 

9 

Aerodynamic  center  location  ahead  of  quarter-chord  at  a 
representative  section  as  ratio  of  average  chord  length 

N.&. 

10 

Quarter-chord  location  ahead  of  elastic  axis  at  a 
representative  section  as  ratio  of  average  chord  length 

n.6. 

11 

Lift  curve  slope 

N.D. 

12 

13 

14 

Drag  coefficient  dg  ~\ 

Drag  coefficient  d^  >  See  Note  #2,  p.  40 

Drag  coefficient  d2 

N.D. 

l$a 

15b 

15c 

15d 

■S  Lag  components  of  four  blade  modes  successively. 

/  Note:  If  less  than  4  blade  modes  are  involved, 

(  the  items  here  must  be  accordingly 

J  limited . 

ft. 

* 

16a 

16b 

16c 

16d 

~n  Flap  components  of  four  blade  modes  in  succession 

C  Note:  As  above.  _ 

ft. 

17a 

17b 

17c 

17d 

\  Rate  of  increase  with  blade  length  of  lag  components 
i  of  blade  modes. 

)  Note:  As  for  15 

N.D. 

18a 

18b 

18c 

18d 

^  Rate  of  increase  with  blade  length  of  flap 
(  components  of  blade  modes. 

)  Note:  As  for  15. 

N.D. 

19a 

19b 

19c 

19d 

N  Twist  components  of  four  blade  modes  in 
)  succession. 

)  Note:  As  for  15. 
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Required  if  left  wing 
is  present,  i.e.  entry 
number  2  on  Card  1  of 
Data  Group  Cl  is  1_. 


The  entry  items 
described  for 
Data  Group  CII 
are  repeated  for 
the  left  wing. 


DATA  GROUP  CIV; 


Required  if  fuselage 
is  present;  .i.e.,  entry 
number  3  on  Card  1  of 
Data  Group  Cl  is  1. 


The  entry  items 
described  for  Data 
Group  C-II  under 
Card  Groups  1 
through  4  are 
repeated  for  the 
fuselage . 


DATA  GROUP  CV: 


Required  if  right  tail 
is  present;  i.e.,  entry 
number  4  on  Card  1  of 
Data  Group  Cl  is  1_. 


Required  if  left  tail 
is  present;  i.e.,  entry 
number  5  on  Card  1  of 
Data  Group  Cl  is  1,. 

Required  if  fin  is  pre¬ 
sent;  i.e.,  entry  num¬ 
ber  on  Card  1  of  Data 
Group  Cl  is  1. 


Entry  items 
described  for 
Data  Group  CII 
are  repeated 
for  right  tail. 

As  above  for  left 
tail . 


As  above  for  the 
f  in . 
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Steady  state  twist  angles 
(including  geometric  twist 
and  collective  pitch) . 
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NOTES: - 


1.  Sign  Convention  for  Rigid  Modes  is: 

Longitudinal  translation  positive  FORWARD 
Lateral  translation  positive  to  the  RIGHT 
Vertical  translation  positive  DOWNWARD 
Roll  translation  positive  RIGHT  SIDE  DOWN 
Pitch  translation  positive  NOSE  UP 
Yaw  translation  positive  RIGHT  SIDE  BACKWARD 

2.  Drag  coefficients  d0,  d,,  d2  are  components 
of  the  Net  Drag  coefficient  (CD)  according 
to  the  formula 

CD  -  d0  +  +  ^ 


where  o<  represents  aerodynamic  angle 
of  attack. 

3.  The  induced  velocity  vector  {vi}  at  a  given 
blade  segment  is  represented  by  a  truncated 
harmonic  series  of  the  azimuth  angle  "V  : 

{vi}  =  {vo]  +  {vlc}co=V'+{vls}SinV+{v2c}c032V+  {v2s}sin2y 


Each  of  the  above  coefficient 
vectors  is  made  up  of  a 
tangential  and  a  radial 
component  each  parallel  to  the 
hub-disc-plane ,  and  a  normal 
component  as  illustrated  in 
the  adjoining  sketch  for  {VOj 
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PROGRAM  OUTPUT 


The  program  generates  a  printed  output  of  several  items  which  fall 
into  three  major  classes  as  described  below. 


1.  Basic  Data:  A  number  of  selected  input  data  and  program  options 

are  printed  out  with  adequate  descriptive  headings 
and/or  other  qualifications  to  help  identify  them. 

2.  Coefficient  Contributions  to  the  coefficient  matrices  by  the 

Matrices:  various  subsystems  and  the  final  matrices  resulting 

therefrom  are  optionally  printed  out  by  rows  with 
row  number  identifications  and  a  descriptive  head¬ 
ing  for  each  complete  matrix.  The  rows  and  columns 
are  associated  with  generalized  forces  and  dis¬ 
placements,  velocities  or  accelerations  correspond¬ 
ing  to  the  system  freedoms  ordered  in  the  following 
sequence: 

(i)  System  rigid  body  freedoms,  if  any; 

(ii)  Airframe  elastic  mode  freedoms,  if  any; 

(iii)  Teetering  freedoms  in  pitch,  yaw  and  roll 
for  the  first  rotor,  as  present; 

(iv)  Blade  mode  freedoms  for  the  first  rotor, 
as  present; 

(v)  Teetering  freedoms  for  the  second  rotor, 
if  any; 

(vi)  Blade  mode  freedoms  for  the  second  rotor, 
as  present;  and 

(vii)  Landing  gear  freedoms,  if  any. 


The  blade  mode  freedoms  are  arranged  with  the 
collective  freedoms  of  all  modes  occurring  first, 
followed  by  the  cyclic  yaw  and  finally  the  cyclic 
pitch  freedoms. 

3.  Problem  Both  the  vibration  and  flutter  solutions  for  a 

Solutions:  given  case  consist  of  a  table  of  complex  roots  of 

the  matrix  equations  and  their  related  physical 
quantities,  and  optionally  the  corresponding 
vectors . 
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The  physical  quantities  related  to  the  roots 
include  frequencies,  a  measure  of  damping  (defined 
differently  for  the  two  solution  types) ,  and  in 
the  case  of  the  flutter  solution,  the  flutter 
velocities . 


Figure  5.  Typical  Landing  Gear  Representation 


SAMPLE  PROBLEM 


INPUT  DATA 


See  following  sheets. 


48 


PROGRAM  OUTPUT 


See  following  sheets. 
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NOTE:  Output  sheets  for  subsequent  cases  have 

been  deliberately  omitted  from  this 
document . 


DISCUSSION  OF  PROGRAM  USAGE 


The  rotor  airloads  calculations  make  use  of  lift  curve  slopes 
input  to  the  program.  It  is  well-known  from  test  data  that 
the  theoretical  value  of  2  it  is  generally  too  high.  A  value  of 
5.73  appears  to  be  a  more  realistic  average  and  is  recommended 
as  a  starting  trial  value  when  no  test  data  for  the  particular 
blade  profile  is  available.  Where  reliable  test  data  exists, 
it  should  form  the  basis  of  the  input  lift  curve  slopes. 

The  induced  airflow  through  the  rotor  disc  does  influence  the 
airloads  calculations,  especially  at  low  advance  ratios. 

Where  no  significant  interference  exists  between  the  airframe 
components  and  the  airflow  at  the  rotor  disc,  it  has  been 
found  adequate  to  include  only  the  airscrew-induced  velocity 
as  a  uniform  axial  velocity  field  calculated  from  a  simple 
momentum  theory. 

The  looping  options  provide  for  changes  to  some  of  the  input 
data  without  affecting  others.  The  concept  involves  some 
basic  assumptions  which  must  be  borne  in  mind  when  the  fea¬ 
tures  are  used.  In  the  rotor  system,  for  example,  the  blade 
modes  do  not  change  when  airspeed  and/or  blade  steady  deflec¬ 
tions  are  varied  under  options  1',  2,  or  3 .  This  is  obviously 
an  approximation  whose  validity  must  be  critically  examined 
in  each  application. 

The  routine  used  in  solving  the  equations  of  motion  involves 
inversion  of  the  stiffness  matrix.  When  rigid  aircraft  trans¬ 
lational  modes  are  included  in  the  analysis  but  no  airframe 
aerodynamics  are  specified,  the  resulting  stiffness  matrix  is 
singular  and  no  solution  can  be  obtained.  In  such  circum¬ 
stances  it  is  advisable  to  specify  the  rigid  translational 
modes  as  airframe  ELASTIC  MODES,  with  the  mode  shapes  at  each 
discrete  point  reflecting  the  actual  motion  under  the  rigid 
mode.  Small  insignificant  modal  frequencies  can  then  be 
specified  for  these  modes  to  keep  the  stiffness  matrix  non¬ 
singular.  Alternatively,  the  landing  gear  facilities  can  be 
used  to  provide  stiffness  terms  in  the  rigid  translational 
modes  without  affecting  the  overall  dynamics  of  the  system. 
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INTRODUCTION 


The  Aeroelastic  Prop/Rotor  Analysis  Program  calculates  the 
following  for  a  single-  or  tandem-rotor  aircraft  in  hover, 
axial  flight,  or  edgewise  flight  (i.e.,  helicopter 
configuration) : 

o  Blade  Airloads  -  Including  downwash,  lift,  drag, 

aerodynamic  pitching  moment,  angle 
of  attack,  and  unsteady  aerodynamic 
parameters,  all  as  a  function  of 
blade  azimuth  and  radial  position 
o  Aerodynamic  Performance  Parameters  -  Including  shaft 

angle  relation  to  free  stream;  ad¬ 
vance  ratio;  rotor  thrust,  drag, 
propulsive  force,  lift,  torque,  and 
horsepower;  and  nondimensionalized 
rotor  thrust,  drag,  propulsive  force, 
lift,  torque,  and  horsepower, 
o  Blade  Response  -  Including  linear  and  angular  deflec¬ 
tions,  velocities  and  accelerations, 
and  shear  and  bending  moment 
distributions 

o  Rotating  and  Fixed  System  Hub  Loads 
o  Pitch  Link  and  Control  System  Actuator  Loads 

The  method  used  employs  an  iterative  technique  between  the  air¬ 
loads  calculations  and  the  blade  response.  Indicators  of  pro¬ 
gram  convergence  are  given  on  the  last  two  pages  of  the  output. 

The  required  input  data  consists  of : 

o  Blade  Physical  Properties 
o  Blade  Airfoil  Characteristics 
o  Trim  Data 

o  Specification  of  Program  Controls 
BLADE  PHYSICAL  PROPERTIES 


The  analysis  idealizes  the  blade  into  a  series  of  bays.  The 
bay  boundaries  are  defined  in  input  locations  290-310.  The 
total  mass  of  each  bay  is  located  midway  between  the  bound¬ 
aries.  The  masses  are  assumed  to  be  connected  by  massless 
beams  which  possess  flap,  lag,  and  torsional  stiffness.  Blade 
property  curves  specify  the  spanwise  distributions  of  blade 
weight,  chordwise  centroid  location,  polar  moment  of  inertia, 
and  flap,  lag,  and  torsional  stiffness  distribution.  The 
lumped  blade  properties  are  obtained  by  inputting  the  distrib¬ 
uted  properties  into  the  D-01  computer  program  and  running 
subroutine  WICK.  A  brief  description  of  D-01,  including  a 
description  of  the  input  and  a  sample  case  based  on  the  dis¬ 
tributed  properties,  is  given  in  a  later  section.  The  output 
of  D-01  subroutine  WICK  may  be  directly  input  into  the  program. 
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BLADE  AIRFOIL  CHARACTERISTICS 


Blade  airfoil  characteristics  for  each  blade  cross  section 
are  defined  by  Cl,  Cp,  and  C^  airfoil  decks.  The  appropri¬ 
ate  deck  must  be  added  to  the  front  of  the  input  data  at  all 
times,  even  when  the  linear  aero  option  is  used.  A  dummy  deck 
should  be  used  since  the  program  reads  the  tables  before  the 
linear  aero  controls.  If  a  blade  with  a  variable  airfoil 
section  is  analyzed,  more  than  one  set  of  airfoil  tables  will 
be  needed  (see  input  description  of  locations  221-223) .  A 
detailed  discussion  of  the  airfoil  tables  is  presented  later. 

Lift  and  moment  y  function  values  associated  with  each  airfoil 
section  are  input  into  locations  145-180.  If  only  a  single 
airfoil  table  is  used,  locations  157  to  180  and  221  to  223 
may  be  set  equal  to  0 . 

TRIM  DATA 

Trim  data  is  usually  obtained  from  either  a  description  of 
wind  tunnel  operating  conditions  or  by  inputting  flight  test 
conditions  into  the  Y-14  trim  program  for  a  single-rotor 
helicopter  or  the  A-97  trim  program  for  a  tandem-rotor 
helicopter.  Program  A-97  is  described  in  a  subsequent  sec¬ 
tion.  Included  is  a  sample  input  and  output  which  serves  as 
the  basis  for  the  sample  C-70  case  presented  later.  The  pro¬ 
gram  input  description  (see  locations  1-21)  illustrates  how 
trim  is  obtained  from  the  A-97  output.  For  a  single-rotor 
helicopter,  locations  15-21  may  be  set  equal  to  0. 

SPECIFICATION  OF  PROGRAM  CONTROLS 

The  computer  program  contains  a  large  number  of  computational 
options;  many  of  these  options  are  used  for  special  purposes 
only.  For  this  reason  a  detailed  discussion  of  all  available 
options  is  not  presented.  The  input  sheets  presented  later 
indicate  usual  values  for  the  program  controls.  Deviation 
from  the  suggested  values  should  be  done  with  due  considera¬ 
tion.  The  other  input  locations  indicate  the  usual  source  of 
the  input  data. 

Location  46  specifies  which  rotor  of  a  tandem  helicopter  is 
to  be  analyzed.  When  analyzing  the  forward  rotor  set  location 
46  equal  to  1.  For  the  aft  rotor  set  location  46  equal  to  2. 
Locations  282  and  283  are  used  to  locate  one  rotor  with  re¬ 
spect  to  the  other  (see  input  description) .  For  a  single¬ 
rotor  helicopter,  set  location  46  equal  to  1  and  locations 
282  and  283  equal  to  0. 

The  most  basic  option  in  this  program  is  a  choice  between 
uniform  downwash  and  nonuniform  downwash.  The  type  of  down- 
wash  used  by  the  program  is  controlled  by  locations  92  and 
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93  (location  93  has  no  effect  in  the  case  of  uniform  down- 
wash)  .  Most  users  will  utilize  either  location  92  equal  to 
0  (nonuniform  downwash)  or  location  92  equal  to  2  (uniform 
downwash)  .  When  location  92  equals  1,  the  program  reads  the 
downwash  as  input  from  locations  586  to  1165.  This  option  is 
useful  for  analyzing  prop/rotors  in  axial  flow,  where  the 
wing  circulation  generates  a  downwash  field  through  the  rotor. 

If  uniform  downwash  is  used,  set  locations  93  to  115  equal  to 
0  EXCEPT  LOCATION  94  WHICH  MUST  BE  SET  EQUAL  TO  1. 

If  the  nonuniform  downwash  option  is  used  then  the  controls 
in  locations  93  to  115  must  be  carefully  specified.  When  a 
tandem-rotor  helicopter  is  being  analyzed,  the  downwash  from 
each  rotor  interferes  with  the  other.  Therefore,  the  down- 
wash  of  both  rotors  must  be  calculated  (deflections  and  loads 
are  calculated  only  for  the  rotor  specified  in  location  46)  . 

If  a  single  rotor  is  being  analyzed,  set  location  93  equal  to 
1  and  only  one  rotor  is  analyzed.  For  a  tandem-rotor  heli¬ 
copter  set  location  93  equal  to  2. 

Locations  94  and  95  are  not  operational  at  the  time  of  this 
writing  and  must  be  set  equal  to  1  and  0,  respectively. 

The  program  calculates  nonuniform  downwash  by  an  iterative 
process;  these  downwash  iterations  are  referred  to  as  loops. 
The  looping  process  is  used  to  insure  compatibility  between 
downwash,  lift,  and  vortex  circulation  strength.  The  looping 
is  performed  separately  for  each  rotor.  Location  96  specifies 
the  number  of  loops  performed  on  the  rotor  being  analyzed  (as 
specified  in  location  46)  and  location  97  specifies  the  number 
of  loops  performed  on  the  interfering  rotor.  Note  that  loca¬ 
tion  97  equals  the  number  of  loops  for  the  interfering  rotor 
and  that  location  96  MINUS  1  is  the  number  of  loops  for  the 
rotor  being  analyzed.  Generally  the  same  number  of  loops  are 
required  for  both  rotors;  hence,  location  96  must  be  1  larger 
than  location  97.  For  a  single-rotor  helicopter  location  97 
is  set  equal  to  0.  When  the  convergence  factors  described  in 
locations  108  to  111  are  used,  most  users  should  be  able  to 
set  location  96  equal  to  1  and  location  97  equal  to  0  and 
achieve  sufficiently  accurate  values  of  circulation  strength. 
Do  not  input  values  larger  than  10  in  locations  96  and  97.  A 
tabulation  of  circulation  strengths  for  each  downwash  loop  and 
a  summary  of  the  thrust  routines  is  provided  to  indicate  the 
convergence  of  the  looping  process. 

The  program  calculates  nonuniform  downwash  by  assuming  that 
each  blade  trails  vortices  in  space.  These  vortices  are 
idealized  into  a  series  of  straight-line  segments.  The  down- 
wash  at  a  given  point  is  calculated  by  applying  the  Biot- 
Savart  law  to  each  straight-line  segment  and  then  summing  for 
all  segments.  Location  98  specifies  the  length  of  vortex  (in 
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rotor  revolutions)  which  is  considered  effective  in  producing 
downwash.  The  input  description  gives  an  equation  by  which 
this  input  value  may  be  calculated. 

Locations  99  and  100  specify  the  inboard  and  outboard  radial 
positions  respectively  from  which  the  vortices  are  trailed. 
Usually  (but  not  necessarily)  the  vortices  are  trailed  from 
the  blade  cutout  and  tip. 

Locations  101  and  102  are  not  operational  at  the  time  of  this 
writing  and  must  be  Set  equal  to  0. 

In  order  to  speed  the  convergence  of  the  downwash  looping 
process,  values  of  circulation  strength  obtained  in  successive 
loops  are  modified:  i.e.,  only  a  portion  of  the  change  is 
added  to  the  value  obtained  in  the  previous  loop.  Locations 
103  and  104  control  the  amount  of  the  change  which  is  added. 
Usually  locations  103  and  104  are  set  equal  to  90.  This  means 
that  10  percent  of  the  change  in  circulation  strength  is  added 
to  the  value  obtained  in  the  previous  loop.  If  a  single  rotor 
is  being  analyzed,  only  location  103  must  be  specified;  for  a 
tandem  both  locations  103  and  104  must  be  specified. 

The  magnitude  of  the  circulation  strength  depends  on  lift. 
Locations  105  and  106  indicate  the  number  of  blade  bays  which 
are  used  to  obtain  the  lift  for  the  outboard  and  inboard 
trailed  vortices  respectively.  Usually  both  locations  105  and 
106  are  set  equal  to  1.  This  means  that  only  the  aerodynamic 
bays  at  the  extremes  of  the  airfoil  (blade  tip  and  cutout)  are 
used  to  determine  lift.  If  both  locations  are  set  equal  to  2, 
then  the  lift  used  will  be  the  average  over  two  bays  and  so 
on.  Locations  105  and  106  do  not  have  to  be  equal. 

The  Biot-Savart  law  predicts  a  1/r2  relationship  between  the 
downwash  at  a  given  point  and  the  distance,  R,  to  a  vortex 
segment.  Thus  if  a  blade  intersected  a  vortex  very  large 
values  of  downwash  would  result.  In  order  to  avoid  this  the 
program  limits  the  downwash  velocity  obtained  in  any  individ¬ 
ual  blade-vortex  interaction  to  0.1  of  the  blade  tipspeed. 

When  location  107  is  set  equal  to  1,  no  upper  limit  is  im¬ 
posed.  The  user  is  advised  to  set  location  107  equal  to  0. 

Experience  has  shown  that  the  circulation  strengths  obtained 
from  the  initial  lift  distribution  are  significantly  higher 
than  the  values  needed  for  lift-circulation  strength-downwash 
compatibility.  Therefore,  in  order  to  speed  convergence  a 
set  of  empirical  factors  has  been  found  which,  when  applied 
to  the  zero-loop  circulation  strengths,  give  the  approximate 
values.  (In  this  way  the  number  of  downwash  loops  can  be  re¬ 
duced.)  These  factors  are  plotted  as  a  function  of  airspeed 
and  are  input  into  locations  108  to  111.  If  zeros  are  input 
into  these  locations,  then  the  factors  are  taken  to  be  0.1. 


102 


Locations  112  to  115  specify  factors  which  modify  the  circu¬ 
lation  strengths  in  all  (including  the  zeroth)  loops.  Usually 
it  is  unnecessary  to  further  modify  the  circulation  strengths 
and  these  factors  should  be  set  equal  to  1.  Do  not  set  loca¬ 
tions  112  to  115  equal  to  0  as  this  will  zero  out  the  down- 
wash. 

The  input  sheets  give  the  source  for  the  required  data  and 
indicate  suggested  values  for  the  program  controls.  Locations 
40,  41,  and  61  show  two  values.  The  values  on  the  left  are  to 
be  used  with  uniform  downwash  (location  92  equal  to  2)  and 
those  on  the  right  should  be  used  with  nonuniform  downwash 
(location  92  equal  to  0)  . 

The  computer  program  uses  an  iterative  technique  between  the 
airloads  and  coupled  flap-pitch  response.  The  number  of 
iterations  performed  is  specified  in  location  48.  It  is  sug¬ 
gested  that  10  iterations  be  used  (the  maximum)  to  provide 
the  greatest  likelihood  of  program  convergence.  Indicators 
of  convergence  are  printed  on  the  last  two  pages  of  output. 

Locations  65  to  69  control  which  parameters  are  included  in 
the  program  printout.  For  most  users  the  values  suggested  on 
the  input  sheet  will  provide  sufficient  information  regarding 
the  intermediate  iterations  (see  sample  output) .  Setting 
location  65  equal  to  1  will  provide  complete  output  for  all 
iterations.  Note  that  lag  and  hub  load  calculations  are  per¬ 
formed  in  the  last  iteration  only. 

The  options  specified  in  locations  586  to  1300  are  not  needed 
by  most  users  and  so  should  be  left  blank;  i.e.,  only  input 
sheets  1  and  2  need  to  be  completed.  In  general,  zeros  need 
not  be  filled  in;  the  program  will  insert  zeros  in  all  blank 
input  locations. 

AIRFOIL  DECKS 

When  large-scale,  high-speed  electronic  computers  became 
available  for  engineering  use,  it  became  possible  to  remove 
many  of  the  rather  dubious  and  extreme  simplifying  assumptions 
which  had  been  necessary  in  the  analysis  of  rotors.  One  of 
the  first  tasks  undertaken  in  making  rotor  calculations  more 
rigorous  was  to  replace  linearized  airfoil  behavior  with  ex¬ 
perimental  data,  including  all  the  awkward  irregularities 
associated  with  stall  and  operation  above  the  critical  Mach 
number. 

This  was  done  by  using  airfoil  decks,  consisting  of  IBM  cards, 
that  defined  an  airfoil's  Cl,  Cq,  and  Cm  characteristics 
versus  angle  of  attack  and  Mach  number.  A  description  of  the 
airfoil  decks  and  a  sample  deck  are  presented  later. 
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Characteristics  of  Airfoil  Tables 


Lift,  drag,  and  pitching  moment  coefficients  for  airfoils  vary 
widely  with  details  of  section  geometry  in  the  range  of  mod¬ 
erate  positive  or  negative  angles  of  attack.  For  angles 
outside  the  range  of  -20  to  +20  degrees  angle  of  attack,  how¬ 
ever,  most  profiles  of  interest  to  the  blade  designer  behave 
very  much  alike.  It  is  therefore  convenient  to  specify  a 
single  function  to  be  used  for  all  different  airfoil  sections 
in  that  range  of  angles.  Within  the  moderate  angle-of-attack 
range,  aerodynamic  coefficients  will  be  specified  in  tabular 
form  against  the  two  arguments:  angle  of  attack  (a)  and  Mach 
number  (M)  ,  which  defines  the  airfoil  deck.  Outside  this 
range  simplified  equations  that  are  written  into  the  rotor 
analysis  are  used  to  define  Cl,  Ce,  and  C^  as  a  function  of 
angle  of  attack. 

The  range  for  which  tabulated  data  is  required  for  lift  co¬ 
efficient  (CL)  is  ±20  degrees;  for  drag  coefficient  (CD)  , 
from  negative  stall  to  positive  stall,  where  the  stall  angles 
are  input  as  aNEG  STALL  and  aSTALL* 

Lift  Coefficient  Deck 


Lift  data  in  the  moderate  angle-of-attack  range  is  specified 
by  a-CL  pairs  at  fixed  Mach  number.  Up  to  10  Mach  numbers 
may  be  used  with  15  C^-a  pairs  per  Mach  number. 

The  following  formulae  and  the  curve  in  Figure  6  define  the 
single  Cl-«  relation  used  for  all  airfoils  in  the  range  not 
covered  by  the  tables  discussed  above.  The  formulae  are  part 
of  the  computer  program,  but  it  will  be  noted  that  four  con¬ 
stants  (K^ ,  K2 ,  K3 ,  K4)  are  also  required.  These  are  normally 
taken  as  one,  but  may  be  adjusted  to  suit  the  needs  of  unusual 
configurations.  Note  that  K2  is  related  to  another  input 
constant  (Cliso)  and  that  they  must  be  mutually  consistent. 
Figure  6  shows  the  Cl_oi  relation  from  +20  to  340  degrees 
(-20  degrees)  for  all  K's  =  1  and  Cliso  = 

The  following  are  equations  of  line  segments  defining  the 
airfoil  lift  coefficient  at  large  angles  of  attack: 

The  value  <5  equals  a  from  20  degrees  through  180  degrees  and 
equals  (360  degrees  -  a)  above  180  degrees  up  to  340  degrees. 
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Figure  6.  Lift  Coefficient  Versus  Angle  of  Attack  for  All 
Airfoils 


(1)  20°<  a <160°  CL  =  (KX  or  K2)  x  (-0.42045519  +  0.08899219 

x  6  -1.5350914  x  10"3  x  <$2  +  7.4254402 
x  10"6  x  63  -8.2868165  x  10"9  x  64) 

(K^  for  20o<a<90°,  K2  for  90°< a xl60° ) 


(2)  160°<  a< 180°  CL  =  (K2)  x  (99.79237-1.2027941  x  6 

+  0.0036020154  x  62)  +  CL180 


(3)  18Q°<a<200° 


(K3)  x  (-99.79237  +  1.2027941 
6  -0.0036020154  <52)  +  CL18Q  ( 


-  160) 
20 


(4)  200°  <ot<  340°  CL  =  (K3  or  K4)  x  (0.42045519  -0.08899219 

x  6  +  1.5350914  x  10“3  x  62  -7.4254402 
x  10-6  x  63  +  8.2868165  x  10-9  x  64) 

(K3  for  200°<a<270° ,  K4  for  270°<a<340°) 


Drag  Coefficient  Deck 

In  the  unstalled  range,  drag  data  are  specified  in  the  form 
CD-M  pairs  for  constant  values  of  a.  Up  to  15  values  of  a 
may  be  used,  each  with  up  to  7  Cd-M  pairs.  In  addition,  a's 
for  negative  and  positive  stall  must  be  specified. 

The  function  applicable  between  the  stall  angles  is  given  in 
the  following  equations  and  the  curve  in  Figure  7. 

The  following  are  equations  of  line  segments  defining  the  air¬ 
foil  drag  coefficient  at  large  angles  of  attack: 

(5)  cD  =  -1.0783242  +  0.10546602  x  v  -0.0022458297  x  v2 

+  2.543888  x  10"5  x  v3  -1.0989072  x  10-7  x  v4 

Q(]°  _ 

where  v  =  a  +  (15°  -  astalx)  (90°  -  astali)  f°X  astall<al90° 
and  v  =  (180°  -a)  for  90°<ct<166° 


106 


(6)  CD  =  5.5012552  -0.060345598  x  6  +  0.00016618434  x  62 

where  6  =  a  for  166<a<_180° 

and  6  =  (360°  -a)  for  180°<a<194° 


(7)  CD  =  -1.0783242  +  0.10546602  x  v  -0.0022458297  x  v2 
+  2.543888  x  10-5  x  v3  -1.0989072  x  10-7  x  v4 


where 
and  v 


v  =  (g  -180°)  for  194°<g<270° 


=  360°  -a  -(345° 
for  270°<_a<aneg 


-aneg 

stall 


stall^ 


(g  -  270° 
aneg  stall 


270° 


) 


Pitching  Moment  Coefficient  Deck 
Moderate  a  Range 

For  -16°<a<16° ,  tabulated  data  must  be  provided,  as 
follows: 


Cm  at  a  =  -16°  (344°),  one  value 

C  at  a  =  16°,  one  value 
m 

Up  to  13  sets  of  Cm-Mach  number  pairs,  corresponding 
to  a  set  of  up  to  13  a's.  Up  to  10  pairs  per  set 
may  be  provided. 

High  a  Range 

The  curve  for  Cm  between  16°<_  x  <_344°  is  presented  in 
Figure  8  and  the  equations  for  the  curve  are  given  below. 
The  first  term  in  the  equation  for  Cm  between  16°  and  30° 
and  between  330°  and  344°  is  zero  for  a  symmetrical  air¬ 
foil  and  not  equal  to  zero  for  cambered  airfcils. 

The  following  are  equations  of  line  segments  defining  the 
airfoil  pitching  moment  coefficient  at  large  angles  of 
attack : 

(8)  le^g^O0  Cm  =  (Cmi6°  +  0.0955)  (3°1~a)  +  (0.22140169 

-0.037698898  x  g  +  0.0014700305  x  g2 
-  2.1980119  x  10"5  x  g3) 
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(9)  30°<  a  <150°  Cm  =  0.16285396  -  0.014395669  x  a  +  1.2695044 

x  10“4  x  a2  -  9.4724852  x  10-7  x  a3 
+  3.7413282  x  10-9  x  a4 

(10)  150° < a_<168°  Cm  =  (-9.8201326  +  0.051520593  x  a 

+  4.2809828  x  10-4  x  a2  -  2.3662806 
x  10-6  x  a3) 

(11)  168°^a<L80°  Cm  =  46.119316  -  0.42395661  x  a  +  0.0011007744 

x  a2  -  6.7548054  x  10-6  x  a3  +  3.2327550 
x  10“8  x  a4 

6  =  360°  -a  and  is  applicable  between  180°  and  344° . 

(12)  180°<a<192°  Cm  =  -46.119316  +  0.42395661  x  6  -0.0011007744 

x  6 2  +  6.7548054  x  10“6  x  6  3  -3.2327550 
x  10-8  x  <54 

(13)  192°<_a<  210°  Cm  =  9.8201326  -  0.051520593  6  -4.2809828 

x  10-4  x  6 2  +  2.3662806  x  10"6  x  63 

(14)  210oj<a<330o  Cm  =  -0.16285396  +  0.014395669  6  -1.2695044 

x  10"4  x  62  +  9.4724852  x  10"7  x  63 
-  3.7413282  x  10-9  x  64 

(15)  330°<_a<344°  Cm  =  (Cm344o  -  0.0955)  (^--^33°)  +  (-0.2214069 

+  0.037698898  x  6  -  0.0014700305  x  62 
+  2.1980119  x  10"5  x  <53 ) 
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A  typical  C_,  C^,  and  airfoil  deck  is  defined  below. 

This  examplS  is  shown  only  to  indicate  how  to  define  the 
decks  and  is  not  meant  to  indicate  actual  airfoil  data. 

C  Deck 

Li 

Title  Card 

180.  180. MAX  POS  ALPHA, MAX  NEG  ALPHA 

0.0  1.0000  1.0000  1.0000  1 . OOOOCONTROL  NOS. 

10  NO. OF  MACH  NUMBERS  FOR  CL  VS  ALPHA 
MACH  NUMBERS 

0.0  0.300  0.400  0.550  0.650  0.750  0.800 

0.850  0.900  1.000 

************************  lift  TABLE  *********************************** 
5  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.  -  0.0 

ALPHA  6 

0.0  6.00  20.000  340.000  360.000 

CL 

-0.010000  0.640000  2.040000  -2.240000  -0.010000 

5  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.300 

ALPHA 

0.0  7.400  20.000  340.000  360.000 

CL 

-0.020000  0.825000  2.146999  -2.219999  -0.020000 

5  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.400 

ALPHA 

0.0  6.500  20.000  340.000  360.000 

CL 

-0.015000  0.755000  2.266999  -2.410000  -0.015000 

5  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.550 

ALPHA 

0.0  3.000  20.000  340.000  360.000 

CL 

-0.025000  0.360000  2.485000  -2.525000  -0.025000 

6  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.650 

ALPHA 

0.0  2.000  4.000  20.000  '  340.000  360.000 

CL 

-0.020000  0.280000  0.555000  2.714999  -3.020000  -0.020000 

5  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.750 
ALPHA 

0.0  1.000  20.000  340.000  360.000 

CL 

-0.020000  0.160000  2.990000  -3.620000  -0.020000 

6  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.800 

ALPHA 

0.0  1.200  20.000  340.000  358.800  360.000 

CL 

-0.045000  0.225000  3.134999  -3.195000  -0.285000  -0.045000 

5  NM.OF  ALPHA-CL  PAIRS  FOR  MACH  NUM.=  0.850 
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ALPHA 

0.0 

3.600 

20.000 

340.000 

360.000 

CL 

-0.045000 

0.795000 

2.106999  -2.540000 

-0.045000 

7  NM.OF 

ALPHA-CL 

PAIRS  FOR  MACH  NUM.= 

0.900 

ALFHA 

0.0 

1.200 

6.200 

20.000 

340.000  358.900 

CL 

-0.015000 

0.075000 

0.800000 

1.766000 

-1.480000  -0.095000 

4  NM.OF 

ALPHA-CL 

PAIRS  FOR  MACH  NUM.= 

1.000 

ALPHA 

0.0 

20 . 000 

340.000 

360 . 000 

CL 

-0.010000 

0.190000 

-1.209999  -0.010000 

360.000 

-0.015000 


C  Deck 
D 


(No  title  card  required. 


attached  to 


Cl  deck) 


POS-NEG  STALL  ANGLES  IN  CD-M  TABLES 
12.000  350.000 

14  NO. OF  ALPHA  VALUES  FOR  CD  VS  M 
ALPHA 


0.0  1.000 

2.000 

4.000 

6.000 

8.000 

12.000  14.000 

350.000 

357.000 

358.000 

359 . 000 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  0. 

0 

MACH 

0.0  0.400 

0.700 

0.810 

0.840 

0 . 877 

CD 

0.011600  0.008600 

0.009200 

0.010100 

0.014900 

0.032800 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  1. 

000 

MACH 

0.0  0.400 

0.750 

0.775 

0.810 

0.850 

CD 

0.011600  0.008800 

0.009300 

0.010300 

0.013900 

0.024400 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  2. 

000 

MACH 

0.0  0.400 

0.625 

0.700 

0.765 

0 . 800 

CD 

0.011600  0.008700 

0.009100 

0.010000 

0.014000 

0.023400 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  4. 

000 

MACH 

0.0  0.300 

0.625 

0.675 

0.725 

0. 800 

CD 

0.011600  0.009200 

0.011100 

0.015300 

0.029000 

0.062900 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  6. 

000 

MACH 

0.0  0.300 

0.400 

0.530 

0.600 

0.650 

CD 

0.043100 

0.011600  0.009200 

0.010550 

0.013500 

0.022100 

7  NUM  OF  H-CD  PAIRS  FOR  ALPHA  =  8.000 


10.000 

360.000 

1.000 

0.105700 

1.000 

0.112700 

1.000 

0.132200 

1.000 

0.202700 

1.000 

0.277700 
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MACH 


0.400 


0.450 


0.0 


0.300 


0.500 


CD 


0.013700  0.011700  0.015800  0.022980  0.034300 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  10.000 

MACH 

0.0  0.200  0.250  0.300  0.350 

CD 

0.018600  0.017400  0.017700  0.023200  0.034400 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  12.000 

MACH 


0.0  0.150  0.200  0.250  0.300 

CD 

0.024600  0.024600  0.027400  0.038700  0.058200 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  14.000 

MACH 

0.0  0.100  0.125  0.150  0.175 

CD 

0.031600  0.031600  0.032100  0.034600  0.038900 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  350.000 

MACH 

0.0  0.250  0.300  0.360  0.410 

CD. 

0.024100  0.022200  0.022400  0.025300  0.030800 

7  NUM  of  M-CD  PAIRS  FOR  ALPHA  =  357.000 

MACH 

0.0  0.500  0.650  0.700  0.750 

CD 

0.012100  0.008800  0.010600  0.016100  0.027000 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  358.000 

MACH 

0.0  0.400  0.700  0.750  0.800 

CD 

0.011600  0.008700  0.010100  0.014500  0.025000 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  359.000 

MACH 

0.0  0.400  0.750  0.775  0.810 


CD 

0.011600  0.008800  0.009300  0.010500  0.015400 

7  NUM  OF  M-CD  PAIRS  FOR  ALPHA  =  360.000 


MACH 

0.0 

0.400 

0.700 

0 . 810 

0.840 

CD 

0.011600 

0.008600 

0.009200 

0.010120 

0.014900 

0.600 

1.000 

0.080100 

0.351200 

0.400 

1.000 

0.052700 

0.427700 

0.400 

1.000 

0.114200 

0.477500 

0.200 

1.000 

0.048900 

0.533700 

0.450 

1.000 

0.040400 

0.259500 

0.825 

1.000 

0.053900 

0.152700 

0.825 

1.000 

0.036400 

0.137700 

0.850 

1.000 

0.025900 

0.117700 

0.877 

1.000 

0.032800 

0.105700 
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C„  DECK 

_ M _ 


CM  ALPHA= 16  DEG  and  344  DEG 
-0.096  0-096 


15  NO. OF  ALPHA  VALUES  FOR  CM  VS  M 
ALPHA 


0.0 

2.000 

4.000  6.000  8.000 

10 . 000 

11.000 

13.000 

16.000 

344.000  349.000  352.000 

356 . 000 

359.000 

360.000 
10  NO. OF 

M-CM  PAIRS 

FOR  ALPHA  =  0.0 

MACH 

0.0 

0.100 

0.200  0.300  0.500 

0 . 600 

0.800 

0.900 

0.950 

1.000 

CM 

0.0 

0.0 

0.0  0.0  0.0  0.0 

0.0 

0.0 

-0.001500  - 

0.002000 

10  NO. OF 

-CM  PAIRS 

FOR  ALPHA  =  2.000 

MACH 

0.0 

0. 300 

0.600  0.700  0.750 

0.800 

0.850 

0.900 

0.950 

1.000 

CM 

-0.002000 

-0.00]  ]00  - 

0.001000  0.000500  0.002200 

-0.001200 

-0.061000 

-0.049000 

-0.035000  - 

0.039500 

10  NO. OF 

M-CM  PAIRS 

FOR  ALPHA  =  4.000 

MACH 

0.0 

0.300 

0.500  0.600  0.650 

0 . 750 

0.800 

0.850 

0 .900 

1.000 

CM 

-0.006000 

-0.004800  - 

0.003000  0.003000  0.004000 

-0.014000 

-0.055500 

-0.072000 

-0.049000  - 

0.012000 

10  NO. OF 

M-CM  PAIRS 

FOR  ALPHA  =  6 . 000 

MACH 

0.0 

0.500  0.600 

0.650 

0.750 

0.800 

0.850 

0.900 

0.925  1.000 

CM 

-0.010000 

-0.008000  0.001000  0. 

001500  -0. 

050000 

-0.070000 

-0.079000 

-0.071000 

-0.046500  -0.080000 

10  NO. OF 

M-CM  PAIRS  FOR  ALPHA  = 

8.000 

MACH 

0.0 

0.300  0.400 

0.500 

0.600 

0. 700 

0.800 

0.900 

0.925  1.000 

CM 

-0.016000 

-0.014000  -0.016500  -0. 

017000  -0. 

012000 

-0.033500 

-0.069000 

-0.077000 

-0.056000  -0.077500 

10  NO. OF  M-CM  PAIRS  FOR  ALPHA  =  10.000 

MACH 

0.0  0.150  0.400  0.500  0.600  0.700  0.800 

0.900  0.925  1.000 
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-0.015000  -0.020000  -0.032000  -0.032000  -0.031000 
-0.078500  -0.057000  -0.074000 

10  NO. OF  M-CM  PAIRS  FOR  ALPHA  =  11.000 

MACH 

0.0  0.150  0.300  0.400  0.600 

0.900  0.925  1.000 

CM 

-0.022500  -0.025500  -0.046500  -0.053000  -0.052500 
-0.092000  -0.073000  -0.086000 

10  NO. OF  M-CM  PAIRS  FOR  ALPHA  =  13.000 

MACH 

0.0  0.150  0.300  0.400  0.600 

0.900  0.925  1.000 

CM 

-0.071000  -0.072000  -0.078000  -0.088000  -0.088000 
-0.112000  -0.099000  -0.105500 

2  NO. OF  M-CM  PAIRS  FOR  ALPHA  =  16.000 

MACH 

0.0  1.000 

CM 

-0.096000  -0.096000 

2  NO.  OF  M-CM  PAIRS  FOR  ALPHA  -  344.000 

MACH 

0.0  1.000 

CM 

0.096000  0.096000 

10  NO. OF  M-CM  PAIRS  FOR  ALPHA  =  349.000 

MACH 

0.0  0.150  0.300  0.400  0.600 

0.900  0.925  1.000 

CM 

0.022500  0.025500  0.046500  0.053000  0.052500 

0.092000  0.073000  0.086000 

10  NO. OF  M-CM  PAIRS  FOR  ALPHA  =  352.000 

MACH 

0.000  0.300  0.400  0.500  0.600 

0.900  0.925  1.000 

CM 

0.016000  0.014000  0.016500  0.017000  0.012000 

0.077000  0.056000  0.077500 

10  NO. OF  M-CM  FAIRS  FOR  ALPHA  =  356.000 


-0.047000  -0.067500 


0.700 


0.800 


-0.066000  -0.081000 


0.700 


0.800 


-0.096500  -0.105500 


0.700 


0.800 


0.066000  0.081000 


0.700 


0.800 


0.033500  0.069000 
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MACH 


0.0 

0.800 

0.500  0.600 

0 . 650 

0.750 

0.800 

0.850 

0.900 

1.000 

CM 

0.006000 

0.004800 

0.003000  -0.003000  -0. 

004000 

0.014000 

0.055300 

0.072000 

0.049000 

0.012000 

10  NO. OF 

M-CM  PAIRS 

FOR  ALPHA  =  359.000 

MACH 

0.0 

0.300 

0.500  0.600 

0. 750 

0.800 

0.850 

0.900 

0.950 

1.000 

CM 

-0.001000 

0.001000  - 

0.001500  -0.003800  0. 

002000 

0.000500 

0.035000 

0.012500 

0.010000 

0 . 019000 

10  NO. OF 

M-CM  PAIRS 

FOR  ALPHA  =  360.000 

MACH 

0.0 

0.100 

0.200  0.300 

0.500 

0.600 

0.800 

0.900 

0.950 

1.000 

CM 

0.0 

0.0 

0.0  0.0  0. 

0 

0.0 

0 . 0 

0.0 

-0.001500  - 

0.002000 
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Figure  11.  Coefficient  of  Pitching  Moment  for  NACA  0011 
(Mod)  Airfoil  Section 


PROGRAM  INPUT 


Dimension  Code 


N-D  Nondimensional 

SL  Slug  -  mass  unit  =  Lb-Sec/Ft 


LOCATION 

1 

2 


3 

4 

5 


6 


SYMBOL 

DESCRIPTION 

DIMENSION 

V 

Flight  path  velocity  (knots)  ,  TAS 

ktB 

a 

Fuselage  water  line  angle  of 
attack  -  measured  from  direction 
of  free  stream  velocity,  not 
horizon  (positive  nose  up)  - 
same  units  and  sign  convention  as 
aero  trim  output. 

deg 

0 

Rotor  speed. 

rpm 

P 

Air  mass  density. 

si/ ft 3 

S 

Speed  of  sound, 

S  =  49.1  V459.6  +  T ’ 

where  T  -  Temperature  in  degrees  F. 

ft/sec 

_ 

Not  used 

--- 

Locations  7  thru  13  define  the  thrust,  collective,  cyclic 
and  flapping  for  the  forward  rotor.  These  locations  are 
also  used  for  an  isolated  rotor  analysis  (as  in  a  single 
rotor  helicopter).  Aft  rotor  data  (Locs.  15-21)  are 
required  when  a  forward  rotor  is  being  analyzed  and  non- 
uniform  rotor  interference  downwash  is  considered  (i.e., 
when  Loc.  46  =  1,  Loc.  92  =  0  or  3  and  Loc.  93  =  3  or  2) . 

7  T  -  Thrust  of  the  forward  rotor  along  lbs 

the  rotor  shaft.  Same  units  and 
sign  convention  as  aero  trim  output. 

Tof  =  TF  in  aero  trim  output,  A97. 


8  0  -  Collective  pitch  angle  for  the  deg 

forward  .rotor.  Equal  to  the  mean 
mechanical  angle  of  attack  at  the 
0.75  r/R  radial  position  -  measured 
from  the  disc  plane  (positive  nose 
up)  .  Same  units  and  sign  convention 
as  aero  trim  program. 
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' 


c-  70 


122 


LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


Note:  Aero  trim  program  A97  prints 
the  collective  angle  at  the  blade's 
center  of  rotation,  THEOF.  There¬ 
fore 

=  THEOF  -  [amount  of  blade  twist 

from  center  of  rotation 
to  .75  r/R] 

or  for  linearly  twisted  blades 

0  -  =  THEOF  -  0.75  0. 
or  t 

where  0  is  the  blade's  linear  twist 

defined  positive  -  nose  down  at  the 
tip. 

Note:  If  thrust  routine  is  used 
(Loc.  49  =  1,  or  2)  program  will  use 
eQ£  as  a  starting  point  for  thrust 

match.  If  thrust  routine  in  not  used 
(Loc.  49  =  0)  the  forward  rotor  collec¬ 
tive  is  set  at  0  f. 


First  harmonic  coefficients  of  control  system  pitch 
angle  for  the  forward  rotor  (Locs.  9  &  10)  -  zero  degree 
azimuth  position  is  trail  aft. 


9 


10 


11 


6lcf  Lateral  cyclic  pitch  angle  for  the  deg 
forward  rotor  -  cosine  coefficient 
of  cyclic  pitch  (positive  -  nose  up 
at  ip  =  0)  .  Same  units  and  opposite 
sign  convention  as  aero  trim  output. 

0lcf  *  -AICF  in  aero  trim  output, 

A97 . 

01  f  Longitudinal  cyclic  pitch  angle  for 
the  forward  rotor  -  sine  coefficient 
of  cyclic  pitch  (positive  -  nose  up, 
at  ij/  =  90)  .  Same  units  and  opposite 
sign  convention  as  aero  trim  output. 

©ls£  =  -BITF  in  aero  trim  output, A97. 


0f  =  0Qf  +  elcf  sin  ip  +  el8f  cos  if, 


8  -  Coning  angle  of  forward  rotor  -  deg 

steady  flapping  angle  (positive  - 
tip  up) .  Same  units  and  sign 
convention  as  aero  trim  program. 
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LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


0Qf  “  AOF  in  aeto  trim  output,  A97. 

First  harmonic  coefficients  of  blade  flapping  angle 
for  the  forward  rotor  (Locs.  12  &  13).  Zero  degree 
azimuth  position  is  trail  aft. 

12  Blcf  Longitudinal  flapping  of  the  for-  deg 

ward  rotor  -  cosine  coefficient  of 
flapping  (positive  -  tip  up  at 
ip  =  0)  .  Same  units  and  opposite 
sign  convention  as  aero  trim  output. 

=  -AIF  in  aero  trim  output,  A97. 


13 


Lateral  flapping  of  the  forward  deg 

rotor  -  sine  coefficient  of  flapping 
(positive  -  tip  up  at  ij;  =  90)  .  Same 
units  and  opposite  sign  convention 
as  aero  trim  output. 


=  -BIF  in  aero  trim  output, 


A97 . 


Bf  =  9of  +  elcf  003  *  +  Blsf  3in  * 
14  -  Not  used. 


Locations  15  thru  21  define  the  thrust,  collective, 
cyclic  and  flapping  for  the  aft  rotor.  Forward  rotor 
data  (Locs.  7  -  13)  are  required  when  an  aft  rotor  is 
being  analyzed  and  non-uniform  rotor  interference 
downwash  is  considered  (i.e.,  when  Loc.  46  =  2, 

Loc.  92  =  0  or  3  and  Loc.  93  =  2  or  3)  . 

15  T  Thrust  of  the  aft  rotor  along  lbs 

the  rotor  shaft.  Same  units  and 
sign  convention  as  aero  trim  output. 

T  _  =  TR  aero  trim  output,  A97. 

Oa 


16  6  Collective  pitch  angle  for  the  deg 

aft  rotor.  Equal  to  the  mean 
mechanical  angle  of  attack  at  the 
0.75  r/R  radial  position  -  measured 
from  the  disc  plane  (positive  nose 
up) .  Same  units  and  sign  convention 
as  aero  trim  program. 


jOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


Note:  Aero  trim  program  A97  prints 
the  collective  angle  at  the  blade's 
center  of  rotation,  THEOR.  There¬ 
fore 

0Qa  =  THEOR  -  [amount  of  blade  twist 

from  center  of  rotation 
to  .75  r/R] 

or  for  linearly  twisted  blades 

0  =  THEOR  -  0.75  0,. 

oa  t 

where  is  the  blade's  linear  twi3t 

defined  positive  -  nose  down  at  the 
tip. 

Note:  If  thrust  routine  is  used 
(Loc.  49  =  1,  or  2)  program  will  use 
0Qa  as  a  starting  point  for  thrust 

match.  If  thrust  routine  is  not  used 
(Loc.  49  =  0)  the  aft  rotor  collec¬ 
tive  is  set  at  0 

oa 

First  harmonic  coefficients  of  control  system  pitch 
angle  for  the  aft  rotor  (Locs.  17  &  18)  -  zero  degree 
azimuth  position  is  trail  aft. 


17  eica  Lateral  cyclic  pitch  angle  for  the  deg 

aft  rotor  -  cosine  coefficient  of 
cyclic  pitch  (positive  -  nose  up 
at  \p  =  0)  .  Same  units  and  opposite 
sign  convention  as  aero  trim  output. 


18 


Isa 


0lca  =  -AICR  in  aero  trim  output, 


A97. 


Longitudinal  cyclic  pitch  angle  for  deg 
the  aft  rotor  -  sine  coefficient  of 
cyclic  pitch  (positive  -  nose  up, 
at  \p  =  90)  .  Same  units  and  opposite 
sign  convention  as  aero  trim  output. 


)^ga  =  -BITR  in  aero  trim  output,  A97 


6a  =  9oa  +  °lca  008  *  +  °lsa  sin  * 
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LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


19  8  Coning  angle  of  aft  rotor  -  deg 

s  steady  flapping  angle  (positive  - 
tip  up) .  Same  units  and  sign 
convention  as  aero  trim  program. 

8  „  =  AOR  in  aero  trim  output,  A97. 
oa 


First  harmonic  coefficients  of  blade  flapping  angle 
for  the  aft  rotor  (Locs.  20  &  21).  Zero  degree 
azimuth  position  is  trail  aft. 


20 


21 


8.  Longitudinal  flapping  of  the  aft  deg 

rotor  -  cosine  coefficient  of 
flapping  (positive  -  tip  up  at 
ij)  =  0)  .  Same  units  and  opposite 
sign  convention  as  aero  trim  output. 

8^ca  =  -AIR  in  aero  trim  output,  A97. 

8,  Lateral  flapping  of  the  aft  rotor  -  deg 

a  sine  coefficient  of  flapping 

(positive  -  tip  up  at  f  =  90)  . 

Same  units  and  opposite  sign 
convention  as  aero  trim  output. 

3lga  =  -BIR  in  aero  trim  output,  A97. 


22  to  39 


e2c  to  Higher  harmonic  coefficients  of 

01OS  control  system  pitch  angle 

(positive  -  nose  up). 


deg 


Higher  harmonic  control  applied  to  forward  rotor  only 
when  Loc.  46  =*  1  and  to  aft  rotor  only  when  Loc.  42  =  2. 

10 

9  *  ^  ^6kc  cos  ^  +  eks  s^n 

k=2 


40 


v.  opt  Control  for  v.  term  in  a  p 
0  =  include  vA  term 
1  =  omit  vft  t£rm 


v^  option  is  used  only  with  non-uniform  downwash 
(Loc.  92  =  0 ,  1  or  3)  . 


41 


CM  opt 


Control  for  CM  option 


0  =  use  Cw  re 


1  =  use 


¥ 


M 

fa^Lref”1"0  ) 


C  +  c 
L  +  LML=0 
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LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


42 


r 


Control  for  the  calculation  of 
zeroth  lift-downwash  iteration 
vortex  circulation  strengths,  r 


=0  use 

r  iji  * 


i=l 


1. 

i 


=1  use  r 

o 


4  7T  CT 

bl  (2-3  m  2)  . 


r 


ic 


r 


Is 


0 

— 6  it  u  CT 
bl(2-3  u  2) 


=2  use  r  =  4 tt  C_ 

o  T 

bl ( 2-3 m  2) 


ric  “  3  a»  co  e 
Sr 

n  -  -  ,  3  a  c 
Is  1  “ 


1b 


4R(2-3  u  2) 


1  c 


-  6ir  M  CT 
bl (2-3 u  2) 


where  r#  «  K(rQ  +  rlc  cos  +  T±s  sin  *) 

b  =  number  of  aerodynamic  blade  elements 
over  which  lift  is  averaged  =  Loc.  105 
or  Loc.  106  (depending  on  whether 

for  tip  or  root) 


L.  =  lift  over  the  ith  aerodynamic  blade 
^  element  at  azimuth  iJj 

U.  =  total  velocity  at  the  ith  station  at 
^  azimuth 

1,  =  length  of  the  ith  aerodynamic  blade 

1  element 


and 
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LOCATION 


44 


SYMBOL  DESCRIPTION  DIMENSION 

k  =  a  constant  factor  =  Loc.  108, 

Loc.  109,  Loc.  110,  or  Loc.  Ill 
depending  on  whether  is  out¬ 
board  or  inboard  or  for  the  for¬ 
ward  or  aft  rotor 

bl  =  number  of  blades 

^lc'^ls  =  1st  harmonic  coefficients  of 
flapping  =  Loc.  12  and  13  for 
forward  rotor  and  Loc.  20  and  21 
for  an  aft  rotor 


H  Control  to  account  for  blade  deflec¬ 

tions  resulting  from  bending-torsion 
coupling 

=  0  Deflections  are  calculated  without 

bending- torsion  coupling  applied 


=  1  or  3  Torsional  deflections  resulting  from 
incremental  changes  in  blade  torsion 
due  to  bending  of  the  torsion  axis  are 
calculated  and  added  to  the  pitch  de¬ 
flections.  These  deflections  are 
proportional  to 


and  hence  are  not  present  for  a 
symmetric  section 


=  2  or  3  Incremental  flap  and  lag  deflections 
resulting  from  bending  of  the  torsion 
axis  and  calculated  and  added  to  the 
blade  deflections 


CpT  a  Control  for  equivalent  linear  aero¬ 
dynamics 


=  0  =  Loc.  51 

bro  =  0 


1 


aro  and  b^  are  calculated  from  airload 
routine 


LOCATION 

SYMBOL 

DESCRIPTION 

DIMENSION 

46 

Rotor 

1  =  analyze  forward  rotor 

2  =  analyze  aft  rotor 

— 

Set  equal  to  1  when  analyzing  an  isolated 
rotor  (as  in  a  single  rotor  helicopter) 

47 

bl 

Number  of  blades  per  rotor 

— 

48 

Number  of  iterations  between  the 
airload  and  coupled  flap  pitch 
routines  -  input  values 
from  0  to  10 

49 

Control  for  thrust  routine 
option. 

0  -  do  not  perform  a  thrust  routine 

1  -  perform  a  thrust  routine.  The 

smallest  value  of  collective 
will  be  found. 

2  -  perform  a  thrust  routine. 

Choose  collective  nearest 
the  initial  value  of  0. 

50 

L.A. 

Control  for  aerodynamic  option 

0  -  use  non-linear,  compressible 

aerodynamics  from  table  look  up 

1  -  use  linear  aerodynamics  where 

CT  =  a  a 

ra^a^2  (  *8) 

CD  -012  +  11  IR/C^T" 

CM  *  0 

51 

a 

00 

Two  dimensional  lift  curve  slope. 

Set  equal  to  5.73 

1/rad 

52 

T.L. 

Tip  loss  -  The  ineffective  length 
of  the  blade  tip,  non-dimensional- 
ized  by  the  blade  radius 

N.D. 

Note:  Be  sure  that  this  length 
does  not  exceed  the  length  of  the 

tip  bay. 
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LOCATION 

53 

54 

55 

56 

57 

58 

59 

60 

61 


SYMBOL  DESCRIPTION 


DIMENSION 


M 


A 


FGC 


F 

G 

K 


AT 


L 


Constant  subtraction  from  to 

reduce  stall  delay  (negative  values 
increase  the  stall  delay) 


Constant  subtraction  from  Aa 

reduce  stall  delay  (negative 
increase  the  stall  delay) 


CM  t0 
values 


Yaw  angle  control 
0  -  calculate  cos  A 
1  -  set  cos  A  =  1 


Controls  selection  of  F&G  functions 
0  -  calculate  F  and  G 
1  -  F  *=  Loc.  57,  G  =  Loc.  58 

Constant  value  of  F  function, 
used  when  Loc.  56  =  1 

Constant  value  of  G  function, 
used  when  Loc.  56  =  1 

Controls  calculation  of  reduced 
frequency 

0,  -  Use  k  =  CQn/2V 

1,  -  Use  k  =  CQ  /-V/e  /2V 

Used  to  limit  Aa  (i.e.,  dynamic 
stall  delay  angle)  in  the 
unsteady  aero  routine 
0,  -  Sets  limit  of  n  radians. 

Any  other  number,  -  Sets  limit 
as  defined  by  the  input  value. 

Control  to  modify_the  slope  used 

to  calculate  the  L  term.  Set 

equal  to 

0,  -  slope  is  n 

1,  -  use  the  equivalent  slope 

i.e.,  • 5 ( Ct  -CT  n) / 

^REF  LO 


(nREF(L) cosA) 


rad 


rad 


N.D. 


N.D. 


deg 
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LOCATION  SYMBOL  DESCRIPTION  DIMENSION 


2,  -  use  the  local  slope 

i-e.,  <cLref+>0i-cLref) 50  ‘ 

3,  -  use  the  local  slope  plus 

yawed  flow,  i .e . , 

^C]^EF+.0$""ClREF*  5°/cosA 

4,  -  slope  is  zero 

62  g  Gravity  acceleration. ■  For  lg  ft 

field  use  32.2  sec2 

63  Steady  induced  velocity  of  ft/sec 

forward  rotor,  positive  down. 

Used  when  Loc.  91  =  0  or  2 

64  VQa  Steady  induced  velocity  of  - 

~  the  aft  rotor,  positive  down. 

Used  when  Loc.  91  =  0  or  1 

65  -  Controls  iteration  output  - 

0,  -  Only  the  final  iteration 
will  be  printed 

1,  -  All  iterations  will  be 

printed 

2,  -  Abbreviated  output.  The 

thrust  routine,  aerodynamic 
parameters  and  pitch  link 
loads  are  printed  out  for 
each  iteration,  with 
complete  printout  for  the 
last  iteration 

Note:  If  Loc.  65  =  1  (all  iterations 

printed)  then  Locs.  66  and  67  will 
affect  output  in  all  iterations. 

If  Loc.  65  =  0  or  2,  then  Loc.  66 
and  67  will  affect  output  for  the 
last  iteration  only. 
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LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


66 


67 

68 

69 


70 

71  to  90 


-  Controls  airload  and  coriolis 

output . 

0,  -  No  output  for  these  loads 

1,  -  Complete  output  (airloads, 

coriolis,  Cl  and  P) 

2,  -  Only  coriolis,  Cl  and  P 

3 ,  -  Only  Cl  and  P 

4,  -  Only  airloads,  Cl  and  P 

(P  —  airloads  perpendicular  to 
the  chord) 

-  Not  used  - 

Set  =  1 

-  Not  used  - 

Set  =  1 

-  Control  for  hub  load  output 

0 ,  -  No  output 

1,  -  Complete  hub  load  output 

-  Not  used  - 

ACm  Pitching  moment  coefficient  - 

correction  applied  at  each 
aerodynamic  blade  element  (from 
the  tip  to  the  cut-out)  to  account 
for  trailing  edge  or  tab  bend. 

This  term  shifts  the  C^  vs  a  curve 
by  the  A  Cm  specified.  A  positive 
ACm  increases  the  nose  up  pitching 
moment . 

Locations  91  to  115  control  the  type  and  magnitude 
of  the  downwash.  When  the  non-uniform  downwash  is 
input,  Loc.  92  =  1,  from  the  table  in  Loc.  586  to 
Loc .  1065,  or  when  uniform  downwash  is  used,  Loc. 
92  =  2,  then  Loc.  93  to  115  have  no  effect  and  may 
be  set  equal  to  zero,  except  Loc.  94  which  must  be 
equal  to  1. 
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LOCATION 

91 


92 


93 


94 


95 


SYMBOL  DESCRIPTION 


Cal  V 
(Tabl§) 


Option  for  the  calculation  of 

uniform  downwash 

0,  -  Obtain  V  c  and  V  from 
.  .  of  oa 

input 

1,  -  Calculate  V  obtain  V 

a  .  of  oa 

from  input 

2,  -  Obtain  V  ^  from  input  and 

calculate  V 

oa 

3,  -  Calculate  V  ,  and  V 

of  oa 


Down  Control  specifying  type  of  in- 
Wash  duced  velocity  to  be  used 

0,  -  Calculate  non-uniform  down- 
wash 

1,  -  Input  downwash  from  a  table 

in  Locs .  586  to  1055 

2,  -  Uniform  downwash  only 

3,  -  Calculate  non-uniform  down- 

wash  and  add  uniform  downwash 


D.W.  Controls  type  of  non-uniform  down- 

wash  used 

1,  -  Self-induced  downwash  only 

(single  rotor) 

2,  -  Total  downwash  (=sum  of  self 

induced  and  rotor  interference) 

3,  -  Rotor  interference  downwash 

only 

NODWLP  Control  to  determine  detail  of 
vortex  structure  for  subsequent 
airload  response  iterations. 

0,  -  Calculate  far  wake  for  0 

iteration  only.  Calculate 
11  filament  mid  and  near 
wake  every  iteration. 

1,  -  Calculate  downwash  for  tip 

and  root  vortex  only,  for 
the  0  iteration. 

2,  -  Calculate  far  wake  for  0 

iterations  only.  Calculate 
11  filament  mid  and  near 
wake  for  0  iterations  only. 

3,  -  Calculate  far  wake  for  every 

iteration.  Calculate  11 
filament  mid  and  near  wake 
for  every  iteration. 

Do  not  use.  Set  equal  to  1. 

NWI  Number  of  near  wake  iteration  loops 

for  the  first  to  tenth  iteration. 


DIMENSION 
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LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


95 (Cont . ) 

96 

97 

98 


99 

100 

101 


Note:  If  equal  to  zero,  the  near 
wake  iteration  loops  may  be  per¬ 
formed  for  the  zero  flap-pitch 
iteration  only. 

Do  not  use.  Set  equal  to  zero. 

LOOPSI  Number  of  lif t-downwash  compatibil¬ 
ity  iterations  for  self-induced  down- 
wash.  Minus  one.  Do  not  exceed  10. 

LOOPRI  Number  of  lif t-downwash  compati-  - 

bility  iterations  for  rotor 
interference  downwash.  Do  not 
exceed  10. 

M  Length  of  trailed  vortex  used  in  N-D 

non-uniform  downwash  calculations, 
expressed  in'  rotor  revolutions. 

M  =  (2  +  Al/R)  /2ir  ti 

where  Al^  =  Loc.  282  for  tandem  rotors 
R 

^1  =  0  for  single  rotor 
R 

Note:  Enter  the  value  of  M  obtained 

from  the  above  equation  raised  to 
the  next  higher  1/10. 

Example:  Calculate  M  =  1.62 

Use  Loc.  98  =  1.7 

Do  not  use  values  of  M  larger  than 

2.0 

a.  Distance  from  center  line  of  hub  N-D 

1  to  inboard  trailed  vortex,  non- 

dimensional  w.r.t.  blade  radius. 

a  Distance  from  center  line  of  hub  N-D 

to  outboard  trailed  vortex,  non- 
dimensional  w.r.t.  blade  radius. 

V  Tangential  and  radial  velocity  - 

option.  If  control  equals 
0,  -  Omit  tangential  and  radial 
velocity  components 
1,  -  Include  tangential  and 
radial  component  of  the 
downwash  in  the  airload 
calculations . 
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LOCATION 

102 


103 


104 


105 

106 

107 

108 


SYMBOL  DESCRIPTION 


a  Not  used 

Set  =  0 


DAMPSI  Damping  factor  applied  to  vortex 

circulation  strengths,  Fiji,  to 
speed  convergence  of  lift-down- 
wash  iterations  for  the  self- 
induced  downwash.  Set  equal  to 
90  (defined  ds  a  percent) . 

DAMPRI  Damping  factor  applied  to  vortex 

circulation  strengths,  Fij;,  to 
speed  convergence  of  lift-down- 
wash  iterations  for  the  rotor 
interference  downwash.  Set 
equal  to  90  (defined  as  a 
percent) . 


Example:  DAMPRI  =  90.  means 
0.1  of  the  change  in  circulation 
strength,  ri|;,  will  be  added  to 
the  previous  F . 


r (n+1) 

after 

damping 


Tn 

after 

damping 


+  0.1  x 


rn+1 

before 

damping 


Tn 

after 

damping 


bT  Number  of  outboard  aerodynamic 

blade  elements  which  the  lift  is 
averaged  to  determine  the  out¬ 
board  vortex  circulation  strength. 


bR  Number  of  inboard  aerodynamic 

blade  elements  over  which  the 
lift  is  averaged  to  determine 
inboard  trailed  vortex  circulation 
strength. 


PRTEST  Control  to  limit  non-uniform  down- 

wash  to  a  maximum  value  of  .IRQ. 

0,  -  Apply  .IRQ  limit  to  downwash 
1,  -  Do  not  apply  limit 


XKINTF  Factor  applied  to  outboard  vortex 

circulation  strength,  Fiji,  in  the 
zeroth  lift  -  downwash  iteration 
of  the  forward  rotor. 


DIMENSION 


pet 


pet 
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LOCATION 

109 

110 

111 

112 

113 

114 

115 

116  to  133 


SYMBOL 


DESCRIPTION 


DIMENSION 


XKINTRF  Factor  applied  to  inboard 

vortex  circulation  strength, 
lip,  in  the  zeroth  lift-down- 
wash  iteration  for  the 
forward  rotor. 

XKINTA  Factor  applied  to  outboard 
vortex  circulation  strength 
I>,  in  the  zeroth  lift-vortex 
downwash  iterations  for  the 
aft  rotor. 

XKINRA  Factor  applied  to  inboard  vortex 
circulation  strength,  Tij>,  in  the 
zeroth  lif t-vortex-downwash 
iteration  for  the  aft  rotor. 


K 


TF 


K 


RF 


K 


TS 


K 


RS 


6kc 

0ks 


Note:  When  location  108  and  111 
are  set  equal  to  0,  the  program 
sets  the  value  at  .1.  If  any 
other  value  is  desired,  simply 
put  the  value  into  the  proper 
location.  See  Figures  12  and  13. 

Multiplication  factor  to  modify 
outboard  vortex  circulation 
strength  for  the  forward  rotor. 

Multiplication  factor  to  modify 
inboard  vortex  circulation 
strength  for  the  forward  rotor. 

Multiplication  factor  to  modify 
outboard  vortex  circulation 
strength  for  the  aft  rotor. 

Multiplication  factor  to  modify 
inboard  vortex  circulation 
strength  for  the  aft  rotor. 

If  non-uniform  downwash  is 
calculated  (Loc.  92  =  0  or  3) 
do  not  set  Locs .  112  to  Loc.  115 
equal  to  zero.  Doing  so  will 
zero  out  the  non-uniform  downwash, 

Second  through  tenth  harmonic 
flap  angle  (positive  -  tip  up) 


Loc . 
Loc . 


116  = 
118  = 


2c 

3c 


Loc. 

Loc. 


117  = 
119  = 


2s 
3s  ' 


deg 


etc , 
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LOCATION 

116  to  133 
(Continued) 


134 

135 

136  to  140 
141 


142 

143,  144 


SYMBOL  DESCRIPTION 


Note:  Use  these  locations  only 
when  Loc.  275  =  8  or  9.  The  steady 
(coning)  and  first  harmonic  are  de¬ 
fined  in  Locs.  11,  12  and  13  for  a 
forward  rotor  and  in  Locs.  19,  20, 

21  for  an  aft  rotor.  The  higher 
harmonic  flapping  coefficients  are 
only  applied  to  the  rotor  being 
analyzed;  i.e.,  forward  rotor  when 
Loc.  46  =  1  and  aft  rotor  when 
Loc.  46  =  2. 

£  Pre-lag  angle 

av  Angle  between  forward  speed  and 

horizontal  plane;  positive  when 
aircraft  is  climbing 

-  Free 

K  Forced  convergence  factor  used  to 

speed  convergence  of  airloads 
response  iterations 

DN  '  dn  +  (dn-i  -  k 

where  =  displacement  used  to 

calculate  N  +  1  iteration 
airloads 

D  =  displacement  based  on  Nth 
iteration  airloads 

°N-1  =  displacements  used  to 

calculate  Nth  iteration 
airloads 

When  K  =  1  forced  convergence  is 
not  employed 

F  Damping  factor  -  the  factor  (1  +  F) 

is  applied  to  the  aero  pitch  damp¬ 
ing  in  the  force  matrix.  This  factor 
is  used  to  speed  convergence. 

Free 


DIMENSION 


deg 

deg 


139 


LOCATION 

SYMBOL 

DESCRIPTION 

DIMENSION 

145 

to 

150 

vL(1) 

Lift  y  function  values  of  first 
airfoil  table  (most  inboard)  for 

Mach  numbers  from  0  to  1. 

N.D. 

151 

to 

156 

(1) 

YM 

Moment  y  function  values  of  first 
airfoil  table  (most  inboard)  for 

Mach  numbers  from  0  to  1. 

N.D. 

157 

to 

162 

(2) 

yl 

Lift  y  function  values  of  second 
airfoil  table  for  Mach  numbers 
from  0  to  1. 

N.D. 

163 

to 

168 

(2) 

ym 

(3) 

Moment  y  function  values  of  second 
airfoil  table  for  Mach  numbers 
from  0  to  1. 

N.D. 

169 

to 

174 

Y  L 

Lift  y  function  values  of  third 
airfoil  table  (most  outboard)  for 
Mach  numbers  from  0  to  1. 

N.D. 

175 

to 

180 

(3) 

ym 

Moment  y  function  values  of  third 
airfoil  table  (most  outboard)  for 
Mach  numbers  from  0  to  1. 

N.D. 

181 

to 

210 

— 

Free 

— 

211 

to 

220 

C-70 

non- 

has  the 
constant 

Free 

capability  of  analyzing  blades  with 
airfoil  sections.  This  is  done  by 

... 

specifying  up  to  3  airfoil  tables  and  the  radial 
positions  at  which  they  apply.  The  program  per¬ 
forms  the  double  table  lookup  for  each  airfoil 
table  and  then  linearly  interpolates  to  obtain 
the  coefficients  for  the  aerodynamic  blade  elements 
in  the  intermediate  positions.  If  a  single  airfoil 
section  is  used  or  if  linear  aerodynamics  is  used 
(Loc.  50  =  l)  set  Loc.  221  to  223  =  0. 


221 


222 

223 


kl 

R 

Radial  position  at  which  most  in¬ 
board  airfoil  table  is  applied. 

Set  equal  to  r/R  of  cutout. 

N.D. 

k2 

R 

Radial  position  at  which  second 
airfoil  table  is  applied. 

N.D. 

k3 

R 

Radial  position  at  which  most  out¬ 
board  airfoil  table  is  applied. 

Set  equal  to  1. 

N.D. 
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LOCATION 

SYMBOL 

DESCRIPTION 

DIMENSION 

Enter  the  airfoil  tables  corres¬ 
ponding  to  these  radial  positions 
in  boxes  at  the  top  of  input  sheet 

No.  1 

224 

— 

Free 

— 

225 

— 

Free 

— 

226  to  265 

— 

Free 

— 

266 

267 

268 

^D2 

AC 

fl^D3 

Increments  of  drag  coefficient 
which  are  added  to  drag  coefficient 
obtained  in  double  table  lookup. 

ACd1  applies  to  the  most  outboard 

table,  ACp2  to  the  second  table,  etc. 

Note :  The  AC^ ' s  are  added  to  the 

drag  coefficients  before  any  inter¬ 
polation  is  performed  for  aerodynamic 
blade  elements  between  the  positions 
where  the  tables  are  applied  (Locs. 

221  to  223) . 

269 

— 

Free 

— 

270 

— 

Free 

— 

271 

R 

Blade  radius 

ft. 

272 

Co 

Reference  chord 

in . 

273 

N .  A . 

Number  of  aerodynamic  blade  elements 
(up  to  20) 

N.D. 

21 A 

N.M. 

Total  number  of  blade  elements  (up 
to  20) 

N.D. 

275 

F.B.B. 

Blade  boundary  control 

— 

=  1 

Rigid  blade 

=  2 

Articulated  blade  -  pitch  housing 
between  flap  and  lag  hinge  -  flap 
hinge  inboard 

=  3 

Articulated  blade  -  pitch  housing 
outboard  of  lag  hinge  -  flap  hinge 
inboard 

141 


LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


276 


a  Distance  from  pitch  arm  to  pitch 

link 


in . 


Positive  toward  center  of  rotation 


277 


Horizontal  hub  offset  (positive  in 
the  direction  of  rotation) . 


in. 


278  b  Pitch  arm,  measured  perpendicular  in. 

to  the  pitch  axis  from  the  pitch 
axis  to  the  pitch  link  (positive 
in  the  direction  of  rotation) . 


142 


Ini 


LOCATION 

279 

Rotor 

280 

281 

282 

283 

284 

285 

286 

287 


SYMBOL 


DESCRIPTION 

Flap-pitch  coupling  angle  - 
angle  through  which  the  "flap" 
hinge  is  rotated. 


aft 
A  l/R 

Ah/R 


K. 


C  =  4  Mp 

•ft  <jj  r 

max 


+  C 


where  C  viscous  torsional  damping 
rate  in  in-lb-sec 
rad 


DIMENSION 

deg 


deg 


Shaft  tilt  of  aft  rotor-positive 
tilted  forward 

Longitudinal  distance  between  rotors, 
non-dimensional  with  respect  to 
blade  radius 

Vertical  distance  between  rotors, 
non-dimensional  with  respect  to 
blade  radius  measured  from  forward 
rotor  to  aft  rotor,  positive  up 

Equivalent  flap  spring 
Equivalent  flap  viscous  damping 

Lag  spring  rate 


Equivalent  lag  torsional  damping 
coefficient 


N.D. 


N.D. 


lb-in/rad 

lb-in-sec 

rad 

lb-in 

rad 

lb-in-sec 

rad 


Mp  pre  load  moment  in  damper  in- lb 

to  lag  natural  frequency  rad/sec 

C  amplitude  of  lagging  motion  radians 
max 


143 


LOCATION 
287  (Cont'd) 


288,  289 
290-315 


316-330 


SYMBOL  DESCRIPTION  DIMENSION 

Note:  The  above  relation  was 
obtained  by  assuming  c  =  Cmax  sin  wt 
and  equating  the  energy  lost  in  one 
cycle  for  an  ideal  torsional  viscous 
damper  (damping  coefficient  =•  C,.)  to 
the  energy  lost  in  1  cycle  for  a 
torsional  viscous  damper  with  a 
constant  coulomb  (pre  load)  damper 
of  Mp  superimposed  and  solving  for 

-  Not  used  - 

xi 

—  Distance  from  rotor  center  to  blade 

element  boundaries  non-dimensionalized 
by  blade  radius,  Loc.  271.  Any  num¬ 
ber  of  boundaries  may  be  used  up  to 
a  maximum  of  21;  i.e.,  maximum  of  20 
elements.  Location  of  boundaries  must 
satisfy  the  following  conditions: 

-Must  have  a  boundary  at  blade  tip; 

i.e. ,  Loc.  290  =  1. 

-Must  have  boundary  at  cutout. 

-Must  have  boundary  at  flap  hinge 
which  is  innermost  boundary.  For 
a  teetering  rotor  flap  hinge  at 
x/R  =  0 . 

-Must  have  boundary  at  lag  hinge 
which  is  adjacent  to  flap  hinge. 

-Number  of  blade  elements  between 
blade  tip  and  cutout  must  equal 
value  in  Loc.  273. 

-Total  number  of  blade  elements 
must  equal  value  in  Loc.  274. 


C  /C 
n-  o 


Blade  chord  non-dimensionalized  by 
reference  chord,  Loc.  272. 


N.D. 


331-335 


Free 


LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


336-360 


361-380 


381-405 


406-420 


M  Lumped  blade  element  masses  for  each 

blade  element  starting  from  the  blade 
tip.  Must  have  a  number  of  entries 
equal  to  value  in  Loc.  274. 

A 6  Incremental  blade  twist  between  masses 

"  1  (masses  are  located  at  centers  of  blade 
elements).  *»tn  -  9t(n+1,  -  8^ 


where  0fc  is  the  built  in  twist  angle 

at  the  nth  mass  measured  from  the 
disc  plane  in  degrees,  positive  nose 
up.  Only  masses  on  the  airfoil 
portion  of  the  blade  are  considered; 
therefore,  n  ranges  from  1  to  (NA-1) 
where  NA  =  Loc.  273. 


The  last  A0 
as 


tn ' 


i .  e .  , 


^  6 1 (NA) 


is  defined 


Y 

on 


t  (NA) 


t  ( -75) 


t  (NA) 


where  0 


and  0, 


are  the 


ut(.75)  t (NA) 

built  in  twist  angles  at  the  .75  r/R 
and  NAth  mass  respectively. 


Distance  from  the  pitch  axis  to  the 
mass  center  of  gravity  (positive 
towards  the  trailing  edge)  for  each 
blade  element,  starting  from  the 
blade  tip. 


Distance  from  mid-chord  to  the 
pitch  axis  (positive-forward)  for 
each  aerodynamic  blade  element, 
starting  from  the  blade  tip. 
Normalized  by  blade  element  chord, 


in . 


N.D. 


C  =  Loc.  316  to  Loc.  335  x  Loc.  272 
n 

421-425  —  Free 


426-450 

451-475 


GJ 

n 


Blade  pitch  inertia  about  the  pitch  lbs-sec- 
axis,  for  each  blade  elemert  in 

starting  from  the  blade  tip. 

.  2 

Torsional  rigidity  between  masses.  lb-m 

When  set  equal  to  0,  the  GJ  is 
assumed  infinite. 


LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


476 

477 

Rotor 

478-481 

482-485 


Rl  Radial  distance  to  the  pitch  link, 

measured  from  the  center  of  rotation. 

Y_  Angular  displacement  of  the  pitch 

"  link  relative  to  the  blade  azimuth 

position,  positive  in  the  direction 
of  rotation. 


R. ,  R^ ,  Radial  distance  to  fixed  system 
R^ ,  R]J  control  actuators,  measured  from 
center  of  rotation. 

Y 1  y o  Azimuth  position  of  the  fixed  system 
Y ^ f  Y4  control  actuators. 

Note:  Loc.  481  and  485  are  used 
only  for  a  4  actuator  control 
system,  with  pivoting  double 
actuators.  If  both  locations  are 
zero,  a  three  actuator  fixed 
control  system  is  assumed. 

When  the  4  actuator  system  is 
analyzed,  1  and  4  are  the  single 
actuators ,  2  and  3  are  the  pivot¬ 
ing  double  actuators . 


i  n . 


deg 


in. 


deg 


LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


486,487  a,  b  Parameters  used  to  expand  the  deg 

region  of  detailed  aerodynamic 
coefficients  defined  in  the  short 
form  CT  airfoil  deck. 

a  -  is  the  low  angle  of  attack 
limit  (tables  currently  stop  at 
20  degrees) . 

b  -  is  the  high  angle  of  attack 
limit  (tables  currently  stop  at 
340  (-20)  degrees) . 

Note:  Detailed  CT  are  currently 

JL 

defined  from  a  =  340  (-20)  to  20 
degrees.  The  airfoil  deck  must 
be  expanded  to  include  the  new 
limits,  but  the  maximum  number  of 
data  points  in  the  table  may  not 
be  exceeded. 

488,  489  c,  d  Parameters  used  to  expand  the  deg 

region  of  detailed  aerodynamic 
coefficients  defined  in  the 
short  form  CM  airfoil  deck. 

c  -  is  the  low  angle  of  attack 
limit  (tables  currently  stop  at 
16  degrees) . 

d  -  is  the  high  angle  of  attack 
limit  (tables  currently  stop  at 
344  (-16)  degrees) . 

Note:  Detailed  CM  are  currently 

defined  from  a  =  344  (-16)  to  16 
degrees.  The  airfoil  deck  must 
be  expanded  to  include  the  new 
limits,  but  the  maximum  number  of 
data  points  in  the  table  may  not 
be  exceeded. 

490  Kz  Control  system  spring  rate.  When  lb/in 

set  equal  to  zero,  K  is  assumed 
infinite . 
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LOCATION 


SYMBOL 


DESCRIPTION 


DIMENSION 


0p  Angle  between  the  pitch  arm  and  the  deg 

disc  plane  for  zero  collective  at 
r/R  =  0.75  -  positive  nose  up 


492  Kg  Stiffness  of  torsional  spring  in. -lb/rad 

between  pitch  housing  and  hub 


493  Cq  Torsional  root  damper  in-lb-sec 

rad 

494  10  Distance  from  root  boundary  to  in 

pitch  bearing 

495  -  Not  used  - 

496-520  E--fn  Flapwise  bending  rigidity  between  lb-in2 

masses.  When  set  equal  to  0,  the 
El  is  assumed  infinite. 

521-540  e  Distance  from  pitch  axis  to  elastic  in 

axis  -  positive  toward  trailing 
edge 

541-565  EIrn  Lagwise  bending  rigidity.  When  set  lb-in2 

equal  to  0,  the  El  is  assumed 
infinite . 
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LOCATION 

SYMBOL 

DESCRIPTION 

DIMENSION 

566-585 

«n 

Distance  from  pitch  axis  to  vertical 
neutral  axis  -  positive  toward 
trailing  edge 

in 

586-1065 

V 

Input  array  for  externally 
calculated  downwash  velocity, 
non-dimensionalized  with  respect 
to  tip  speed  ) positive  down) . 

Use  only  when  location  92  =  1. 

N.D. 

1066-1085 

— 

Free 

— 

1086-1110 

Amn 

Cabled  mass . 

lb-sec^/ 

in 
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SIGN  CONVENTIONS 


a.  Sign  Convention  for  Rotor  Hub  Separation 
(Aft  Rotor  shown  as  First  Rotor) 

.Tof 


Disc  Plane 


Hater  Line 


Relative  Hind 


Figure  14.  Fuselage  Geometry 


Blade 


44  > 


Outboard  Vortex 


Inboard 

Vortex 


cb 


£ 


Figure  15.  Vortex  Geometry 


Figure  19.  Radial  Force  Component  From  Flap  Foreshortening 
Sign  Convention 
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Response  Conventions :  The  sign  conventions  for  the 
forces  and  displacements  in  coupled  flap-pitch  and 
coupled  lag  are  shown  below: 


Flap  Pin 


V  Lag  Pin 


<!/'  / 


»«7 


Figure  21.  Coupled  Flap-Pitch  Sign  Convention 


Flap 

Fin  , 


(JV 


:  The  convention 
appllea  to  either 
rotor.  Direct Iona 
remaining  conalatent 
with  the  direction 
of  rotation 


# 

Figure  22.  Uncoupled  Lag  Sign  Convention 


e.  Hub  Load  Sign  Convention 


Figure  23.  Rotating  and  Fixed-System  Hub  Loads  Sign  Convention 
(Shaft  Axis  System) 


*«A 


Figure  24.  Steady  Loads  Relative  to  Airstream 
(Wind  Axis  System) 


DESCRIPTION  OF  OUTPUT  DATA 


The  program  output  is  available  in  several  versions  depending 
on  problem  requirements.  See  locations  65  to  69.  Described 
below  is  the  output  data  in  each  of  the  program  sections. 

I  Input  Data  -  The  input  location  number,  n,  and  the  value 
are  printed  out. 

ANALYSIS  OUTPUT 

SYMBOL  NAME  DESCRIPTION  DIMENSION 

Loc  n  Item  n  Input  location  number  n 

Value  Value  input  in  location  n  See  input 

description 

II  Initial  Deflection  -  No  printout. 

III  Downwash  -  The  uniform  downwash  velocity  normalized  by 
the  tipspeed  for  both  rotors  is  printed.  This  is  done 
for  all  downwash  options. 

If  the  program  internally  calculates  nonuniform  downwash 
(Location  92  =  0  or  3) ,  then  the  following  is  printed: 

1.  Tip  and  root  circulation  strength  normalized  by  rotor 
angular  velocity  versus  blade  station  and  azimuth 
position  for  one  or  both  rotors,  depending  on 
location  93. 

2.  Downwash  velocity  normalized  by  tipspeed  versus  blade 
station  and  azimuth  position  for  one  or  both  rotors, 
depending  on  location  93. 

3.  Tip  and  root  circulation  strengths  normalized  by  rotor 
angular  velocity  for  the  last  downwash  loop,  before 
and  after  damping,  versus  azimuth  position  for  one  or 
both  rotors,  depending  on  location  93. 

4.  The  total  downwash  velocity  normalized  by  the  tip- 
speed  versus  blade  station  and  azimuth  position  for 
the  rotor  being  analyzed. 

Total  downwash  =  self-induced  downwash 

+  rotor  interference  downwash 

5.  A  summary  of  the  thrust  routines  used  in  the  downwash 
calculations  with  complete  details  for  thrust  routine 
used  in  the  zeroth  and  last  downwash  loops. 

If  nonuniform  downwash  is  specified  in  locations  586  to 
1065;  i.e.,  location  92  =  1,  then  the  following  is 
printed : 
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1.  Downwash  velocity  normalized  by  tipspeea  versus  blade 
station  and  azimuth  for  the  rotor  being  analyzed. 

Following  the  downwash  routine  output,  the  program  prints 
the  results  of  the  airloads  and  coupled  flap-pitch 
routines  for  each  iteration.  The  specific  airload  param¬ 
eters  which  are  printed  are  controlled  by  the  values  in 
locations  65-67.  Use  of  the  thrust  routine  is  controlled 
by  locations  49  and  275. 

IV  Airloads  -  Complete  output  of  the  airloads  routine 
includes : 


1.  Thrust  routine. 

2.  Unsteady  aerodynamics  parameters  at  each  blade  sta¬ 
tion  and  azimuth  position. 

3.  Airload  perpendicular  to  the  blade  chord  at  each 
blade  station  and  azimuth. 

4.  Vibratory  aerodynamic  pitching  moment  per  unit  length 
about  the  pitch  axis. 

5.  Classical  aerodynamic  parameters  (lift,  drag, angles 
of  attack,  etc.)  at  each  blade  station  and  azimuth 
position . 

6.  Normal  and  radial  components  of  thrust  in  harmonic 
form  at  each  blade  station. 

7.  Component  of  drag  in  the  disc  plane  and  the  total 
aerodynamic  pitching  moment  in  harmonic  form  for  each 
blade  station. 

8.  Aerodynamic  coefficients,  Cl,  Cq,  C^,  as  obtained 
from  the  double  table  look-up  routine  (after  inter¬ 
polation)  at  each  blade  station  and  azimuth. 

THRUST  ROUTINE 


ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION 


i 


STA  Blade  station  number 


DIMENSION 


ATn  1st  DEL  T  Steady  thrust  along  the  rotor  lb 

shaft  at  blade  station  n  for  the 
first  value  of  collective  used 
by  the  thrust  routine . 

Positive  up 
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ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 

ATn  Last  DEL  T  Steady  thrust  along  the  rotor  lb 

shaft  at  blade  station  n  for  the 
last  value  of  collective  used  by 
the  thrust  routine. 

Positive  up 

04  THETA  AT  blade  mechanical  angle  of  radians 

STA.  4  attack  at  station  4 

Positive  leading  edge  up 

Tq  THRUST  Total  rotor  thrust  along  rotor  lb 

shaft 

Positive  up 
UNSTEADY  AERODYNAMICS 

x  PSI  Blade  azimuth  -  positive  in  deg 

direction  of  blade  rotation 

F  F  Value  of  in-phase  Theodorsen  ND 

deficiency  function 

G  G  Value  of  Theodorsen  deficiency  ND 

function  phased  90°  behind  F 

aEQUIV  ALPHA  EQUIV  Dynamic  angle  of  attack  -  radians 

angle  of  attack  including 
Theodorsen  deficiency  func¬ 
tions  F  and  G 

a  ref  l  ALPHA  Quasi-static  angle  of  attack  radians 

REF (L)  adjusted  for  dynamic  lift 

stall  delay 

“REF  M  ALPHA  Quasi-static  angle  of  attack  radians 

REF (M)  adjusted  for  dynamic  moment 

stall  delay 

CLref  CL (REF)  Quasi-static  lift  coefficient  - 

based  on  oiref  l 

CLo_i8o  CL(0-180)  Static  lift  coefficient  at  zero  - 

angle  of  attack,  if  in  reverse- 
flow  region  angle  of  attack  at 
180° 

L  L  BAR  Lift  due  to  blade  pitching  lb 

velocity 
Positive  up 
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ANALYSIS 

OUTPUT 

SYMBOL 

SYMBOL 

DESCRIPTION 

DIMENSION 

*  aBE 

ALPHA  DOT 

Time  derivative  of  ALPHA  be 

rad 

Positive  leading  edge  up 

sec 

• 

H 

H  DOT 

Velocity  of  air  normal  to  disc 

17 

plane 

sec 

Positive  up 

Note:  If  uniform  downwash  only  is  used;  i.e.,  location 
92  =  2,  the  downwash  table  versus  blade  station  and 
azimuth  is  filed  with  the  value  of  the  uniform  downwash. 


< 

o 

t-h 

VZF 

Uniform  downwash  velocity  for 
forward  rotor 

Positive  down 

ND  WRT 
tipspeed 

Vo 

s 

VZS 

Uniform  downwash  velocity  for 
aft  rotor 

Positive  down 

ND  WRT 
tipspeed 

X 

PSI 

Blade  azimuth  position 

Positive  in  direction  of 
rotation 

deg 

LOOP 

Number  of  lif t-downwash 
compatibility  iterations 

— 

r 

GAMMA 

Vortex  circulation  strength 
BEFORE  refers  to  before  damping 
AFTER  refers  to  after  damping 
Positive  counterclockwise 
looking  in  direction  of  blade 
velocity 

ND  WRT 
blade 
angular 
velocity 

i 

STA 

Blade  station  number 

V 

VA 

Downwash  velocity 

Positive  down 

ND  WRT 
tipspeed 

04 

THETA 

Blade  mechanical  angle  of 
attack  at  station  4 

Positive  leading  edge  up 

radians  - 

To 

THRUST 

Total  rotor  thrust 

Positive  up 

lb 

NO  ITER 

Number  of  attempts  required 

— 
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ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

ATn 

DELT 

Steady  thrust  at  blade  station 
n  along  rotor  shaft.  FIRST 

DELT  designates  conditions  for 
the  first  value  of  collective 
used  by  the  thrust  routine. 

LAST  DELT  designates  conditions 
for  last  value  of  collection 
used  by  thrust  routine. 

Positive  up 

lb 

SL 

L-B 

One-half  of  the  lift  slope  for 

the  L  term.  CL 

da 

1 

rad 

COS 

COS (YAW) 

Cosine  of  yaw  angle 

— 

PHI 

Dynamic  induced  blade  angle 
of  attack 

Positive  leading  edge  up 

d  CL 
da  a=0 

REF  SLOPE 

Lift  slope  for  a  equal  to  zero 

1 

rad 

AIRLOAD  PERPENDICULAR  TO  CHORD 

AL 

“  “  “ 

Total  airload  perpendicular  to 
blade  chord 

Positive  up 

lb 

r/R 

RAD 

Radial  position  of  mass  ND  WKT  blad< 

radians 

AL0 

STDY 

Steady  value  of  total  airload 
perpendicular  to  blade  chord 
Positive  up 

1 

BAYL 

Length  of  aerodynamic  bay 

in. 

AL-AL0 

1 

“ 

Vibratory  airload  perpendicular 
to  blade  chord  per  unit  length 
Positive  up 

lb/in . 

AMPL 

AMPL 

Amplitude  of  airload  perpen¬ 
dicular  to  chord  per  unit 
length  -  equal  f o  half  peak- 
to-peak  value 

lb/in . 
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ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

AL 

1  O 

STDY 

Steady  value  of  total  airload 
perpendicular  to  chord  line 
per  unit  length 

Positive  up 

lb/in . 

Q 

Vibratory  aerodynamic  pitching 
moment  about  the  pitch  axis 
per  unit  length 

Positive  leading  edge  up 

in. -lb 
in. 

AMPL 

AMPL 

Amplitude  of  the  aerodynamic 
pitching  moment  about  pitch 
axis  per  unit  length 

in . -lb 
in . 

Qo/i 

STDY 

Steady  value  of  aerodynamic 
pitching  moment  about  the 
pitch  axis  per  unit  length 
Positive  leading  edge  up 

in . -lb 
in . 

CLASSICAL  AERODYNAMICS 

■‘•n 

STATION 

Blade  station 

— 

r/R 

x/R 

Radial  position  for  station  n 

ND  to  blade 
radius 

X 

PSI 

Blade  azimuth  position 

Positive  in  direction  of 
blade  rotation 

deg 

0 

THETA 

Blade  mechanical  angle  of 
attack 

Positive  leading  edge  up 

rad 

a 

ALPHA 

Blade  total  angle  of  attack 
Positive  leading  edge  up 

rad 

PHI 

Induced  blade  angle  of  attack 
Positive  leading  edge  up 

rad 

AP 

LAMBP 

Velocity  component  of  the  air 
normal  to  disc  plane 

Positive  up 

ND  to 
tip'-peed 

At 

LAMBT 

Velocity  component  of  the  air 
in  the  disc  plane 

Positive  toward  trailing  edge 

ND  to 
tipspeed 

MACH 

MACH 

Mach  number  based  on  U 

— 
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ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

U 

U 

Blade  resultant  velocity 

ND  to 
tipspeed 

dfx 

DFX 

Bay  radial  force  due  to 
thrust 

Positive  away  from  hub 

lb 

dT 

T 

Bay  thrust  along  rotor  shaft 
Positive  up 

lb 

dH 

H 

Bay  horizontal  force 

Positive  toward  trailing  edge 

lb 

dQ 

Q 

Total  aerodynamic  pitching 
moment  about  pitch  axis 

Positive  nose  up 

in. -lb 

dL 

L 

Bay  lift  perpendicular  to 
relative  wind 

Positive  up 

lb 

dD 

D 

Bay  drag  force  parallel  to 
relative  wind 

Positive  in  direction  of  wind 

lb 

dM 

M 

Aerodynamic  pitching  moment 
about  aerodynamic  center 
Positive  leading  edge  up 

in. -lb 

NORMAL  AND  RADIAL  COMPONENTS  OF  THRUST  IN  HARMONIC  FORM 

K 

K 

Harmonic  number 

fxn 

FXN 

Total  aerodynamic  radial  force 
at  nth  bay 

Positive  away  from  hub 

lb 

dT 

DEL  T 

Bay  thrust  perpendicular  to 
disc  plane 

Positive  up 

lb 

dH 

DEL  H 

Component  of  the  aerodynamic 
force  in  the  disc  plane  per¬ 
pendicular  to  the  pitch  axis 
Positive  toward  trailing  edge 

lb 

dQ 

DEL  Q 

Aerodynamic  pitching  moment 
about  the  pitch  axis 

Positive  leading  edge  up 

in . -lb 

ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 


CL 


CD 


CM 


RESULTANT  Resultant  (of  harmonic 

THRUST  analysis)  of  bay  thrust 


RESULTANT  Resultant  (of  harmonic 

DRAG  analysis)  of  component  of 

aerodynamic  force  in  disc 
plane  perpendicular  to  pitch 
axis 

Positive  toward  trailing  edge 


CL  Dynamic  lift  coefficient 

obtained  by  multiplying 
dynamic  lift  slope,  d  CL 

bY  aEQUIV  d“ 


CD  Coefficient  of  drag  obtained 

from  double  table  look-up 
including  Mach  number  and 

aREF  M 


CM  Dynamic  pitching  moment  co¬ 

efficient  obtained  by  using 
“REF  M  an^  accounting  for  the 
shift  in  center  of  pressure 
due  to  stall 


lb 


ND 


ND 


ND 


V  Aerodynamic  Performance  Parameters  -  Following  the  air¬ 
loads  certain  aerodynamic  performance  parameters  are 
printed  out.  Listed  below  are  the  aerodynamic  parameters 
and  their  corresponding  symbols. 


ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

aR 

ALPHA  R 

Shaft  angle  relative  to  stream 
velocity 

deg 

V 

VKNOTS 

Freestream  velocity 

knots 

V 

Rfi 

V/R*OMEGA 

Nondimensional  freestream 
velocity 

- 

U 

MU 

Advance  ratio 

- 

“270 

ALPHA 

AT  270 

Angle  of  attack  at  4)  =  270° 

To 

TS 

Rotor  thrust  perpendicular  to 
disc  plane 

lb 
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ANALYSIS 

OUTPUT 

SYMBOL 

SYMBOL 

DESCRIPTION 

DIMENSION 

Ts 

TAU  S 

Steady  torque 

lb-ft 

RHP 

RHP 

Rotor  horsepower 

hp 

H 

H 

Drag  in  disc  plane 

lb 

X 

X 

Propulsive  force 

lb 

L 

L 

Lift  force 

lb 

X 

X 

Ratio  of  propulsive  force 

- 

L 

L 

to  lift  force 

L 

L 

Lift  force  nondimensional 

— 

De 

DE 

to  equivalent  drag 

L 

L 

Nondimensionalized  lift 

— 

qd2a 

QD2SIG 

X 

X 

Nondimensionalized  propulsive 

- 

qd2q 

QD2SIG 

force 

P 

P 

Nondimensionalized  rotor 

qd^a 

QD2SIGV 

horsepower 

Si 

CT 

Nondimensionalized  steady 

a 

SIG 

thrust 

Cip 

CTP 

Nondimensionalized  lift 

a 

SIG 

CH 

CH 

Nondimensionalized  drag 

a 

SIG 

Cp 

CP 

Nondimensionalized  rotor 

a 

SIG 

horsepower 

Ax 

X 

Radial  foreshortening  due  to 
flap  and  lag  motion 

Positive  toward  tip 

In  i-x  form  for  aerodynamic 
stations 

in. 

AC 


C 


Coriolis  inertia  force 
Positive  toward  leading  edge 


lb 


ANALYSIS 

OUTPUT 

SYMBOL 

SYMBOL 

DESCRIPTION 

DIMENSION 

Ax 

XK 

Radial  foreshortening  due  to 
flap  and  lag  motion 

Positive  toward  tip 

In  harmonic  form  for  aero¬ 
dynamic  stations  only 

in. 

VI  After  the  airloads,  the  program  performs  the  response 

calculations.  The  coupled  flap-lag-pitch  routine  results 
include  blade  linear  and  angular  deflections  in  both  the 
undeflected  blade  and  disc  plane  coordinate  systems.  It 
also  produces  the  blade  shears  and  moments  in  both  the 
undeflected  blade  and  blade  coordinate  systems. 

COUPLED  FLAP-PITCH  RESPONSE  in  harmonic  form  in  unde¬ 
flected  blade  coordinate  system  (cosine  coefficients  on 
top) 


ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 


vz 


My 


E 


Z 


MX 


e 


vz 


MY 


BETA 


Z 


MX 


THETA 


Vertical  shear  in  Z  direction  lb 

in  undeflected  blade  coordinate 
system 
Positive  up 


Flap  bending  moment  about  Y 
axis  in  undeflected  blade 
coordinate  system 
Positive  blade  tip  up 

Blade  flap  angle  in  unde¬ 
flected  blade  coordinate 
system 

Positive  blade  tip  up 

Blade  flap  deflection  in  un¬ 
deflected  blade  coordinate 
system 
Positive  up 

Torsional  moment  about  X 
axis  in  undeflected  blade 
coordinate  system 
Positive  leading  edge  up 

Blade  pitch  angle  in  unde¬ 
flected  blade  coordinate 
system 

Positive  leading  edge  up 


in . -lb 


radians 


in. 


in. -lb 


radians 
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ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

Vy 

VY 

Chordwise  shear  in  Y  direction 
in  undeflected  blade  coordinate 
system 

Positive  toward  leading  edge 

lb 

MZ 

MZ 

Lag  bending  moment  about  Z 
axis  in  undeflected  blade 
coordinate  system 

Positive  trailing  edge  in 
tension 

lb 

C 

ZETA 

Blade  lag  angle  in  undeflected 
blade  coordinate  system 

Positive  blade  leading 

radians 

Y 

Y 

Blade  lag  deflection  in 
undeflected  blade  coordinate 
system 

Positive  toward  leading  edge 

in. 

Note:  For  a  blade  which  is  articulated  in  lag,  the  lag 

boundary  is  at  the  lag  hinge  and  so  one  less  mass  station 
is  considered. 


COUPLED  FLAP-LAG-PITCH  RESPONSE  as  a  function  of  blade 
station  and  azimuth  in  the  undeflected  blade  coordinate 
system  including  amplitudes  and  steadies.  The  transfor¬ 
mation  from  harmonic  to  time  history  form  is  done  in  the 
last  iteration  only. 


ANALYSIS  OUTPUT 
SYMBOL  SYMBOL 


Vz  VZ 


My  MY 


f  BETA 


DESCRIPTION  DIMENSION 

Vibratory  vertical  shear  in  lb 

Z  direction  in  undeflected 
blade  coordinate  system 
Positive  up 

Vibratory  flap  bending  moment  in. -lb 

about  Y  axis  in  undeflected 
blade  coordinate  system 
Positive  blade  tip  up 


Vibratory  flap  angle  in  unde-  radians 

fleeted  blade  coordinate 

system 

Positive  tip  up 
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ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 

Z  Z  Vibratory  flap  deflection  in  in. 

undeflected  blade  coordinate 
system 

Positive  tip  up 

Mx  MX  Vibratory  torsional  moment  in. -lb 

about  X  axis  in  undeflected 
blade  coordinate  system 
Positive  leading  edge  up 

0  THETA  Vibratory  pitch  angle  in  un-  radians 

deflected  blade  coordinate 
system 

Positive  leading  edge  up 

Vy  VY  Vibratory  chordwise  shear  in  lb 

Y  direction  in  undeflected 
blade  coordinate  system 
Positive  toward  leading  edge 

Mz  MZ  Vibratory  lag  bending  moment  in. -lb 

about  Z  axis  in  undeflected 
blade  coordinate  system 
Positive  trailing  edge  in 
tension 

5  ZETA  Vibratory  lag  angle  in  unde-  radians 

fleeted  blade  coordinate 
system 

Positive  blade  leading 

Y  Y  Vibratory  lag  deflection  in  the  in. 

undeflected  blade  coordinate 
system 

Positive  toward  trailing  edge 

Vx  VX  Vibratory  radial  shear  in  X  lb 

direction  in  undeflected  blade 
coordinate  system 
Positive  away  from  hub 

A 0  DELTA  THETA  Incremental  torsional  deflection  rad 

due  to  bending  of  the  torsional 
axis  in  undeflected  blade 
coordinate  system 
Positive  leading  edge  up 


mmmm 


ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 

AY  DELTA  Y  Incremental  lag  deflection  due  in. 

to  bending  of  the  torsional 
axis  in  undeflected  blade 
coordinate  system 
Positive  toward  trailing  edge 

AZ  DELTA  Z  Incremental  flap  deflection  due  in. 

to  bending  of  the  torsional  axis 
in  undeflected  blade  coordinate 
system 
Positive  up 

COUPLED  FLAP-LAG-PITCH  RESPONSE  as  a  function  of  blade 
station  and  azimuth  in  the  blade  coordinate  system  in¬ 
cluding  amplitudes  and  steadies.  Only  stations  outboard 
of  the  lag  hinge  are  considered.  These  results  are 
calculated  in  the  last  iteration  only. 


ANALYSIS  OUTPUT 
SYMBOL  SYMBOL 


DESCRIPTION 


DIMENSION 


Vibratory  radial  shear  in  X  lb 

direction  in  blade  coordinate 
system 

Positive  away  from  hub 

Vibratory  chordwise  shear  in  lb 

Y  direction  in  blade  coordinate 
system 

Positive  toward  leading  edge 

Vibratory  vertical  shear  in  lb 

Z  direction  in  blade  coordinate 
system 
Positive  up 

Vibratory  torsional  moment  in. -lb 

about  X  axis  in  blade 
coordinate  system 
Positive  leading  edge  up 

Vibratory  flap  bending  moment  in. -lb 

about  Y  axis  in  blade 
coordinate  system 
Positive  leading  edge  up 

Vibratory  lag  bending  moment  in. -lb 

about  Z  axis  in  blade 
coordinate  system 

Positive  trailing  edge  in  tension 


i 


COUPLED  FLAP -LAG-PITCH  linear  and  angular  deflections, 
velocities,  and  accelerations  in  disc  plane  coordinate 
system  as  a  function  of  blade  station  and  azimuth. 


ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

Y 

Y 

Lag  deflection  in  disc  plane 
coordinate  system 

Positive  toward  trailing  edge 

in . 

Z 

Z 

Flap  deflection  in  disc  plane 
coordinate  system 

Positive  up 

in . 

0 

THETA 

Torsional  angular  deflection 
in  disc  plane  coordinate 
system 

Positive  leading  edge  up 

rad 

e 

BETA 

Flap  angular  deflection  in 
disc  plane  coordinate  system 
Positive  tip  up 

rad 

s 

ZETA 

Lag  angular  deflection  in 
disc  plane  coordinate  system 
Positive  blade  leading 

rad 

VII  After 

the  response  for  the  last  iteration  is 

printed  the 

program  calculates  the  hub  and  control  loads .  The  calcu¬ 
lated  quantities  include  hub  loads  due  to  one  blade 
expressed  in  the  rotating  disc  plane  coordinate  systems 
and  the  total  hub  loads  due  to  all  blades  expressed  in 
both  the  fixed  and  rotating  disc  plane  coordinate  sys¬ 
tems.  Also  calculated  are  the  pitch  link  and  actuator 
loads . 


HUB 

LOADS  FOR 

ONE  BLADE  ONLY 

ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

k 

K 

Harmonic  number 

*xb 

X  XB 

Root  radial  force  due  to 
flap  foreshortening  in 
rotating  disc  plane 
coordinate  system 

Positive  away  from  hub 

lb 

168 


ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

Xyb 

X  YB 

Root  radial  force  due  to 
coriolis  acceleration  in 
rotating  disc  plane 
coordinate  system 

Positive  away  from  hub 

lb 

fxb 

F  XB 

Root  radial  force  due  to  thrust 
in  rotating  disc  plane  coordi¬ 
nate  system 

Positive  away  from  hub 

lb 

Vxb 

V  XB 

Total  root  radial  force  (XXB 
+  XYB  +  FXB)  in  rotating  disc 
plane  coordinate  system 
(steady  value  also  contains 
ax) 

Positive  away  from  hub 

lb 

Fx 

FX 

Radial  hub  load  for  one  blade 
in  rotating  dj  sc  plane  coordi¬ 
nate  system 

Positive  away  from  hub 

lb 

fy 

FY 

Tangential  hub  load  for  one 
blade  in  rotating  disc  plane 
coordinate  system 

Positive  toward  leading  edge 

lb 

Fz 

FZ 

Axial  hub  load  for  one  blade 
in  rotating  disc  plane  coordi¬ 
nate  system,  in  direction  of 
rotor  shaft  axis 

Positive  in  the  direction  of 
positive  thrust 

lb 

Mx 

MX 

Pitching  moment  at  center  of 
rotation  for  one  blade  in 
rotating  disc  plane  coordinate 
system 

Positive  leading  edge  up 

in. -lb 

My 

MY 

Flapping  moment  of  center  of 
rotation  for  one  blade  in 
rotating  disc  plane  coordinate 
system 

Positive  blade  tip  up 

in . -lb 
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ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 

Mz  MZ  Torque  at  center  of  rotation  '  in. -lb 

for  one  blade  in  rotating 
disc  plane  coordinate  system 
(same  as  fixed  disc  plane 
system  torque) 

Positive  trailing  edge  in  torsion 
TOTAL  ROTATING  HUB  LOADS  AT  ROTOR  CENTER  FOR  ALL  BLADES 
ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 

Fxr  FXR  X  component  of  load  trans-  lb 

mitted  to  hub  by  all  blades 
in  rotating  disc  plane  coordi¬ 
nate  system  (X  axis  rotates 
with  blade  being  analyzed) 

Positive  away  from  hub 

Fyr  FYR  Y  component  of  load  trans-  lb 

mitted  to  hub  by  all  blades 
in  rotating  disc  plane  coordi¬ 
nate  system 

Positive  toward  leading  edge 

Fzr  FZR  Axial  component  of  load  trans-  lb 

mitted  to  hub  by  all  blades  in 
rotating  disc  plane  coordinate 
system 

Positive  in  direction  of 
positive  thrust 

Mxr  MXR  Moment  about  X  axis  trans-  in. -lb 

mitted  to  hub  by  all  blades 
in  rotating  disc  plane  coordi¬ 
nate  system  (X  axis  rotates 
with  blade  being  analyzed) 

Positive  leading  edge  up 

Myr  MYR  Moment  about  Y  axis  trans-  in. -lb 

mitted  to  hub  by  all  blades 
in  rotating  disc  plane  coordi¬ 
nate  system 
Positive  blade  tip  up 

Mzr  MZR  Torque  at  center  of  rotation  in. -lb 

due  to  all  blades  in  rotating 

disc  plane  coordinate  system 
(same  as  fixed  dis^  plane  system) 

Positive  trailing  edge  in  tension 
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ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

Fxf 

FXF 

X  component  of  load  trans¬ 
mitted  to  hub  by  all  blades  in 
fixed  disc  plane  coordinate 
system 

Positive  in  direction  zero 
degree  azimuth  position 

lb 

Fyf 

FYF 

Y  component  of  load  trans¬ 
mitted  to  hub  by  all  blades 
in  fixed  disc  plane  coordi¬ 
nate  system 

Positive  in  direction  of 
90-degree  azimuth  position 

lb 

Fzf 

FZF 

Axial  component  of  load  trans¬ 
mitted  to  hub  by  all  blades  in 
fixed  disc  plane  coordinate 
system 

Positive  in  direction  of 
positive  thrust 

lb 

Mxf 

MXF 

Moment  due  to  all  blades 
about  Y  axis  in  fixed  disc 
plane  coordinate  system 

Positive  using  right-hand  rule 

in . -lb 

Myf 

MYF 

Moment  due  to  all  blades 
about  Y  axis  in  fixed  disc 
plane  coordinate  system 

Positive  using  right-hand  rule 

in . -lb 

MZf 

MZF 

Torque  at  center  of  rotation 
due  to  all  blades  about  Z 
axis  in  fixed  disc  plane  system 
Positive  using  right-hand  rule 

in. -lb 

Fx 

FXF  A 

Steady  force  in  X  direction  of 
wind  axis  coordinate  system 
(X  axis  is  downwind) 

Positive  in  direction  of 
positive  X 

lb 

Fy 

FYF  A 

Steady  force  in  Y  direction  of 
wind  axis  coordinate  system 
(same  as  force  in  Y  direction 
of  fixed  disc  plane  coordinate 
system) 

Positive  in  direction  of 
positive  Y 

lb 
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ANALYSIS  OUTPUT 

SYMBOL  SYMBOL  DESCRIPTION  DIMENSION 


Fz  FZF  A  Steady  force  in  Z  direction  lb 

of  wind  axis  coordinate 
system 

Positive  in  direction  of 
positive  Z 

Mx  MXF  A  Steady  moment  about  X  axis  in. -lb 

in  wind  axis  coordinate 
system 

Positive  using  right-hand  rule 

My  MYF  A  Steady  moment  about  Y  axis  in. -lb 

in  wind  axis  coordinate 
system  (same  as  moment  about 
Y  axis  in  fixed  disc  plane 
coordinate  system) 

Positive  using  right-hand  rule 

Mz  MZF  A  Steady  moment  about  Z  axis  in. -lb 

in  wind  axis  coordinate 
system 

Positive  using  right-hand  rule 
RHP  RHP  Rotor  horsepower 

PITCH  LINK  AND  ACTUATOR  LOADS 


DIMENSION 

lb 

lb 

lb 

lb 


ANALYSIS 

SYMBOL 

OUTPUT 

SYMBOL 

DESCRIPTION 

DIMENSION 

PSI 

Blade  azimuth  position 

deg 

P 

F 

Vibratory  pitch  link  load  in 
time  history  form 

Positive  in  compression 

lb 

P1 

PI 

Vibratory  load  in  first  actu¬ 
ator  in  time  history  form 
Positive  in  compression 

lb 

p2 

P2 

Vibratory  load  in  second 
actuator  in  time,  history 
form 

Positive  in  compression 

lb 

p3 

P3 

Vibratory  load  in  third 
actuator  in  time  history 
form 

Positive  in  compression 

lb 

Note:  Maximum  positive,  maximum  negative,  and  the 
amplitude  (=  1/2  MAX  POS  -  MAX  NEG)  are  also  indicated 
for  the  time  histories  of  the  vibratory  pitch  link  and 

actuator  loads . 

PERFORMANCE  PARAMETERS  BASED  ON  HUB  LOADS 


B'ollowing  the  hub  loads  calculation,  the  program  reeval¬ 
uates  some  of  the  performance  parameters  in  terms  of  the 
hub  loads . 

VIII  PITCH  LINK  LOAD  SUMMARY  PAGE 

For  all  output  options  a  summary  of  the  pitch  link  loads 
is  provided.  The  vibratory  pitch  link  loads  are  printed 
in  time  history  form  in  15-degree  azimuth  increments  for 
all  iterations.  The  corresponding  amplitudes  and  steady 
values  are  given.  In  addition,  the  last  iteration  vibra¬ 
tory  pitch  link  loads  in  time  history  form  are  given  in 
5-degree  azimuth  increments. 

IX  CONVERGENCE  SUMMARY  PAGE 

For  all  output  options  a  summary  of  the  program  conver¬ 
gence  is  given  on  the  last  page  of  the  output.  The 
quantities  described  are  given  in  the  table  below  and  are 
listed  for  each  iteration. 
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ANALYSIS 

SYMBOL 

e4 

x/R 

To 
«4  0 

“4  90 

a4  180 

a4  270 

Z 

e 


OUTPUT 

SYMBOL 

COLLECTIVE 

X/R 

THRUST 
PSI  =  0 

PSI  =  90 

PSI  = 180 


PSI  = 270 


Z  TIP  AMPL 


THETA 
TIP  AMPL 


ROOT  FBM 
AMPL 


DESCRIPTION 


Blade  mechanical  angle  of 
attack  at  station  4 

Radial  position  at  which 
64  is  specified 


Total  rotor  thrust 

Blade  angle  of  attack  at 
station  4  and  0-degree 
azimuth  position 
Positive  leading  edge  up 

Blade  angle  of  attack  at 
station  4  and  90-degree 
azimuth  position 
Positive  leading  edge  up 

Blade  angle  of  attack  at 
station  4  and  180-degree 
azimuth  position 
Positive  leading  edge  up 

Blade  angle  of  attack  at 
station  4  and  270-degree 
aximuth  position 
Positive  leading  edge  up 

Blade  tip  flapping  deflection 

amplitude 

Always  positive 

Blade  tip  torsional  deflec¬ 
tion  amplitude 
Always  positive 

Root  flap  bending  moment 

amplitude 

Always  positive 


DIMENSION 

deg 

ND  with 
respect  to 
radius 

lb 

deg 

deg 

deg 

deg 

in. 

deg 

in . -lb 
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SAMPLE  PROGRAM 
INPUT  AND  OUTPUT 
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ROTOR  ANALYSIS  pTFTi 
PROGRAM  C^O  BO-105 
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Similar  output  for  Z.  THETA,  BETA,  ZETA 
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COMPREHENSIVE  ROTOR  ANALYSIS 
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SAMPLE  BLADE  PROPERTY  IDEALIZATION  PROGRAM 
INPUT  AND  OUTPUT 


The  Blade  Property  Idealization  Program  converts  the  contin¬ 
uous  blade  property  distributions  into  a  system  of  discrete 
point  masses  connected  by  massless  elastic  elements.  The 
blade  properties  considered  are  blade  mass,  chordwise  centroid 
polar  mass  moment  of  inertia,  flapwise  flexural  stiffness, 
chordwise  flexural  stiffness,  and  torsional  rigidity. 

A  sample  input  and  output  are  presented.  This  output  was  used 
to  obtain  the  input  for  the  sample  C-70  program.  Also  pre¬ 
sented  are  the  distributed  blade  property  curves . (Figures  25 
through  32)  which  provided  the  input  data. 

Complete  sheets  1,  1A,  2,  3,  4,  and  5. 

1.  Basic  Control  Sheet  -  Locs  1  through  30. 

LOC  DEF  DIM.  DESCRIPTION 


01 

Prog 

Not  used. 

Set  equal  to  1 

02 

Flap 

Bdry 

- 

Not  used. 

Set  equal  to  1 

03 

Lag 

Bdry 

- 

Not  used. 

Set  equal  to  1 

04 

R 

ft 

Blade  radius 

05 

Co 

in. 

Blade  constant  chord 

06 

eo 

in. 

Not  used. 

Set  equal  to  0 

07 

b 

in. 

Not  used. 

Set  equal  to  1 

08 

KZ 

lb/in. 

Not  used. 

Set  equal  to  0 

09 

Ke 

lb-in. 

rad 

Not  used. 

Set  equal  to  0 

10 

n 

rpm 

Not  used. 

Set  equal  to  1 

11 

0)  _ 

S 

HZ 

Not  used. 

Set  equal  to  0 
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9 


m 


-m 


LOC 

DEF 

DIM. 

DESCRIPTION 

12 

“f 

Hz 

Not  used. 

Set  equal  to  1 

13 

A  co 

Hz 

Not  used. 

Set  equal  to  1 

14 

n 

- 

Number  of  mass  stations 

15 

Inp 

- 

Set  equal  to  1 

16 

D.O.F. 

- 

Not  used. 

Set  equal  to  1 

17 

nB 

- 

Total  number  of  blade  masses 

18 

E 

- 

Not  used. 

Set  equal  to  0 

19 

K5 

lb-in. /rad 

Not  used. 

Set  equal  to  0 

20 

e6 

in . 

Flap  hinge  offset 

21 

e5 

in. 

Lag  hinge  offset 

22 

CS 

lb-sec 

in. 

Not  used. 

Set  equal  to  0 

23 

IF 

chug  in.2 

Not  used. 

Set  equal  to  1 

24 

chug  in.2 

Not  used. 

Set  equal  to  1 

25 

BL 

- 

Not  used. 

Set  equal  to  1 

26 

M 

He 

chugs 

Not  used. 

Set  equal  to  0 

27 

Mhl 

chugs 

Not  used. 

Set  equal  to  0 

28 

Ke 

in .  lb 
rad 

Not  used. 

Set  equal  to  0 

29 

C6 

chugs/sec 

Not  used.  . 

Set  equal  to  0 

29 


2. 

Subroutine  WICK  Input  Description  -  Locs  100-1000 

LOC 

DEF 

DIM. 

DESCRIPTION 

100 

Wick 

- 

Control  for  use  of  subroutine 
WICK.  Set  equal  to  0 

101 

M 

chugs 

Mass  at  bay  centers 

Set  equal  to  1 

102 

M 

chugs 

Cabled  mass 

Set  equal  to  0 

103 

I 

chug-in. 2 

Pitch  inertia 

Set  equal  to  1 

104 

Y 

in. 

Distance  from  pitch  axis  to 
bay  mass  center 

Set  equal  to  1 

105 

Blank 

- 

Not  used. 

Set  equal  to  0 

106 

ei6 

lb-in.2 

Flapwise  flexural  rigidity 

Set  equal  to  1 

107 

EI5 

lb  .-in. 2 

Lagwise  flexural  rigidity 

Set  equal  to  1 

108 

GJ 

lb-in. 

Torsional  rigidity  of  bay 

Set  equal  to  1 

109 

EA 

lb 

Elastic  axis  rigidity 

Set  equal  to  0 

110 

- 

- 

Use  for  numbering  cases 

Set  equal  to  1 

111 

R 

in. 

Blade  radius 

112-161 

r  or  r/R 

in.  or  N.D. 

Blade  boundary  locations 

NOTE :  These  blade  boundaries 

must  correspond  to  those  used 
in  C-70  Locs  290-315.  Therefore, 
last  boundary  is  flap  hinge 
location  and  must  be  in  inches; 
i.e.,  same  as  Locs  20,  171,  371, 
471,  581,  681,  and  901.  The  next  to 
the  last  boundary  is  lag  hinge 
location  and  must  be  in  inches; 
i.e.,  same  as  Locs  21  and  791. 
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LOC 

DEF 

DIM. 

162 

in. 

170 

- 

- 

171 

e6 

in. 

172-269 

r  or 

r/R 

in.  or  N.D. 

271-369 

Cl) 

lb/in. 

370 

- 

- 

371 

e6 

in. 

372-419 

r  or 

r/R 

in .  or  N . D . 

421-469 

0) 

lb/in. 

470 

- 

- 

471 

e(3 

in . 

472-519 

r  or 

r/R 

in.  or  N.D. 

521-569 

Ip  X 

102 

lb-sec2-in . 

571 


572 

Co 

in. 

573 

- 

in. 

581 

e8 

in. 

582-629 

r  or  r/R 

in.  or 

631-679 

7 

in. 

DESCRIPTION 


Lag  hinge  location 

Use  for  numbering  cases 
Set  equal  to  1 

Flap  hinge  location 

Blade  boundary  locations 

Running  weight 

Cabled  weight  case 
Use  for  numbering  cases 
Set  equal  to  1 

Flap  hinge  location 

Set  equal  to  zero  if  not  using 

cabled  weight 

Blade  boundary  locations 

Running  cabled  weight 

Pitch  inertia  case 
Using  for  numbering  cases 
Set  equal  to  1 

Flap  hinge  location 

Blade  boundary  locations 

Running  pitch  inertia 

Chordwise  centroid  case 
Use  for  numbering  cases 
Set  equal  to  1 

Blade  constant  chord 

Pitch  axis  offset  (aft  of  L.E.) 

Flap  hinge  location 

Blade  boundary  locations 

Mass  offset  from  leading  edge 

Flap  El  case 

Use  for  numbering  cases 

Set  equal  to  1 


680 


LOC  DEF 

DIM. 

DESCRIPTION 

681  eg 

in. 

Flap  hinge  location 

682-730  r  or  r/R 

in.  or  N.D. 

Blade  boundary  locations 

731-780  Elg 

X  1CT6 

lb-in.2 

Flap  rigidity  (running) 

790 

Lag  El  case 

Use  for  numbering  cases 

Set  equal  to  1 

791  e? 

in. 

Lag  hinge  location 

792-840  r  or  r/R 

in.  or  N.D. 

Blade  boundary  locations 

841-890  EI? 

X  10"6 

lb-in. 2 

Lag  rigidity  (running) 

900 

— 

Torsional  rigidity  case 

Use  for  numbering  cases 

Set  equal  to  1 

901  e3 

in. 

Flap  hinge  location 

902-950  r  or  r/R 

in.  or  N.D. 

Blade  boundary  locations 

951-  GJ 

1000  X  10-6 

lb-in. 2 

Torsional  rigidity  (running) 

NOTE:  Input  9999999.  in  the  location  immediately 

following  the  last  entry  in  each  group  of  input 
data.  See  sample  input  sheet.  Distributed 
properties  are  input  beginning  at  the  inboard 
end  of  the  blade. 


Figure  26.  Effective  Centroidal  Axis  Location  of  CH-47  Rotor  Blade 


Figure  27.  Effective  Polar  Mass  Moment  of  Inertia  About  Pitch 
Axis  of  CH-47  Rotor  Blade 


BLADE  STATION,  r/R  R  =  360  IN 


220 


221 


Figure  30.  Torsional  Stiffness  Distribution  of  CH-47  Rotor  Blade 


Re 

quired  Distribul 

:ion 

Control 

0 

100 

Blade  Weight. 

m 

1 

101 

1.  1 

Cabled  Weight. 

2 

102 

Pitch  Inertia 

i 

3 

103 

1  l. 

Chordwise  Centroid,  y 

4 

104 

j  i. 

5 

105 

0. 

Flap  El 

EI « 

6 

106 

i. 

Lag  El 

El  o 

7 

107 

i  1. 

Torsion  GJ. 

GJ 

8 

108 

1.  1 

Elastic  Axis. 

EA 

9 

109 

_1  L_..J 

5 

.16 

6 

.17 

7 

.18 

8 

.19 

9 

.20 

lr(in)  or  r/R 
360. 
0.99 
0.95 


8.0 

9999999. 


20.21 

21.22 

22.23 

23.24 


Boundar 


b26 .27 
27.28 
b28 . 29 
b29 . 30 


b30.31 
b31 .32 
b32 . 33 
b33 . 34 
b34 . 35 


b35 

.36 

b36 

.37 

b37 

.38 

b38 

.39 

b39 

.40 

b40 

.41 

b41 

.42 

CM 

43 

.43 

b43 

.44 

b44 

.45 

b45 

.46 

b46 

.47 

b47 

.48 

b48 

.49 

b49 

.50 

b50 

.51 

Lag  Hinge  162 


r(in)  or  r/ 


5 


or  N.D. 


m* 


Ufr.-LLl'^ 


SHEET  5 


Loc 


Lag  El  Case  No.  I  790 


Loc  No. 


Torsional  Rigidity  Case  I  900 


14  804 

15  805 


r(in)or 

r/R 


29.5 

0.125 

0.125 

0.187 

0.267 


0.295 

0.342 

0.381 

0.519 

0.537 


0.672 

0.672 

0.761 

0.761 

0.960 


0.9 

1.0 

999999J  858 


21 

22 

23 


24  814 

25  815 


26 

27 

28 

29 

30 


31  821 

32  822 

33  823 


34  824 

35  825 


826 
827 

38  828 

39  829 


Eli] 


[•HU 


44  834 

45  835 


EIxlO! 
lb  in2 


365 

365 


1603. 

1603. 

19999999 


r(in)or 

r/R 


8. 

0.125 
0.165  953 

0.1875  954 

0.294  955 


.339 
.342 

908  0.381 

909  0.519 


GJxlO'6 
Loc  lb  in^ 


202 
202 
953  I  190 


95 


82.1 


Pill] 


(■Tail 


999999 


21 

921 

22 

922 

23 

923 

24 

924 

tUiNil 

mi 


28 

29 

30 

928 

929 

930 

31 

931 

32 

932 

33 

933 

34 

934 

35 

935 

937 


38  938 

39  939 


tVl 


PrT'l 


44 

944 

45 

vim 

46 

946 

47 

947 

48 

948 

49 

949 

50 

950 

VVI 


9999999 


Km 


H-TU 


p*‘*l 


hmt  r~i 
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NO. 

SYM 

UNIT 

DEFINITION 

27 

CF 

ft 

Forward  rotor  blade  chord  length  - 
obtained  from  the  projection  of  the 
leading-edge  and  trailing-edge  span 
lines  to  the  center  of  rotation 

28 

bF 

- 

Number  of  blades  in  the  front  rotor 

29 

eF 

ft 

Front  rotor  blade  flapping  hinge 
offset  (usually  written  as  egF) 

30 

°1F 

Increment  in  solidity  at  the  tip 
(forward  rotor)  based  on  a  linear 
variation  from  the  reference  root 
(input  27,  Cp)  to  the  tip: 
o ip  =  bp  (Crpjp  -  Cp)/icRp 

31 

0TWF 

deg 

Angle  of  total  twist  of  forward 
rotor  blade  from  reference  root  at 
the  centerline  of  rotation  to  the 
tip 

Positive  leading  edge  up  at  blade 
tip 

32 

XCF 

- 

Forward  rotor  radius  ratio  for 
blade  cutout 

33 

k6F 

“ 

Rate  of  change  of  front  rotor  blade 
pitch  angle  with  rotor  blade  flap 
angle.  Equal  to  -  tan  63 

34 

iF 

deg 

Angle  of  incidence  of  the  front 
rotor  shaft  in  the  helicopter  X  -  Z 
plane 

Positive  into  the  wind 

35 

VTF 

ft/sec 

Front  rotor  blade  tipspeed 

36 

IF 

slugs-ft^ 

Mass  moment  of  inertia  of  the 
forward  rotor  blade  about  the  flap¬ 
ping  hinge  (usually  written  as  IBF) 
Obtain  from  D-01 

37 

^WF 

ft-lb 

Weight  moment  of  the  forward  rotor 
blade  about  the  flapping  hinge  in  a 

horizontal  position  (gMg) .  Obtain 
from  D-01 
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NO. 

SYM 

UNIT 

DEFINITION 

38 

bitf 

deg 

Front  rotor  longitudinal  cyclic 
pitch  angle  due  to  trim  devices 
other  than  longitudinal  control. 
Generally  a  function  of  dynamic 
pressure 

44 

1F 

ft 

Distance  from  the  helicopter  eg  to 
the  projection  of  the  front  rotor 
hub  on  the  X  axis 

45 

hF 

ft 

Distance  from  the  helicopter  eg  to 
the  front  rotor  hub  measured 
parallel  to  the  helicopter  Z  axis 

46 

dF 

ft 

Lateral  distance  from  helicopter 
eg  to  the  front  rotor  shaft 

47 

rr 

ft 

Rear  rotor  blade  radius  measured 
from  the  center  of  rotation 

48 

CR 

ft 

Rear  rotor  blade  chord  length  - 
obtained  from  the  projection  of  the 
leading-edge  and  trai ling-edge  span 
lines  to  the  center  of  rotation 

49 

bR 

- 

Number  of  blades  in  the  rear  rotor 

50 

eR 

ft 

Rear  rotor  blade  flapping  hinge 
offset  (usually  written  as  eQ_) 

pK 

51 

°1R 

Increment  in  solidity  at  the  tip 
(rear  rotor)  based  on  a  linear 
variation  from  the  reference  root 
(input  48,  CR) 

52 

0TWR 

deg 

Angle  of  total  twist  of  rear  rotor 
blade  from  reference  root  at  the 
centerline  of  rotation  to  the  tip. 
Positive  leading  edge  up  at  blade 
tip . 

53 

XCR 

- 

Rear  rotor  radius  ratio  for  blade 
cutout 

54 

k6R 

- 

Rate  of  change  of  rear  rotor  blade 
pitch  angle  with  rotor  blade  flap 

angle.  Equal  to  -  tan  53 
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NO.  SYM  UNIT 


117  a± 

to  to 

125  o9 


126 

to 

134 


eTWi 

to 

eTW9 


DEFINITION 

The  solidity  (based  on  local  chord 
bc^/uR)  value  corresponding  to  the 
above  Xn  stations.  Leave  blank  if 
solidity  is  input  in  equation  form. 

The  twist  value  (STW  )  corresponding 
to  the  specified  Xn  location.  Leave 
blank  if  twist  is  input  in  equation 
form. 


NOTE:  C-70  input  descriptions  indir-te  how  to  calculate  the 
required  trim  data  and  identify  the  appropriate  A-97 
output  needed. 
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■  1  3 .  ABSTRACT 
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on  program  usage  are  presented. 

The  methods  contained  in  this  report  are  intended  to  be  used  by  designers  to 
calculate  with  improved  accuracy,  the  dynamic  and  aeroelastic  response  characteris¬ 
tics  of  rotor  powered  V/STOL  aircraft.  The  essential  new  feature  of  these  methods 
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calculations  are  essential  if  a  high  level  of  confidence  is  to  be  had  in  the  results. 
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